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An  Algorithm  and  Associated  Computer  Subroutine  For  Calculating  Flow  Through 
a Horizontal  Ceiling/Floor  Vent  in  a Zone-Type  Compartment  Fire  Model 


Leonard  Y,  Cooper 


ABSTRACT 


An  algorithm  and  associated  computer  subroutine  is  presented  for  calculating 
the  effects  on  two -layer  compartment  fire  environments  of  the  quasi -steady 
flow  through  a horizontal  vent  connecting  two  spaces.  The  two  spaces  can  be 
either  two  inside  rooms  of  a multi-room  facility  or  one  inside  room  and  the 
outside  ambient  environment  local  to  the  vent.  The  description  of  the  flow 
through  the  vent  is  determined  by  combining  considerations  of  1)  the 
unidirectional - type  of  flow  driven  by  a cross-vent  pressure  difference  and, 
when  appropriate,  2)  the  exchange  - type  of  flow  induced  when  the  fluid 
configuration  across  the  vent  is  unstable,  i.e.,  when  a relatively  cool,  dense 

gas  in  the  upper  space  overlays  a less  dense  gas  in  the  lower  space.  In  the 

algorithm,  calculation  of  the  rates  of  flow  exchange  between  the  two  spaces  is 
based  on  a previously  reported  theory.  Characteristics  of  the  geometry  and 
the  instantaneous  environments  of  the  two  spaces  are  assumed  to  be  known  and 
specified  as  inputs.  The  outputs  calculated  by  the  algorithra/subroutine  are 
the  rates  and  the  properties  of  the  vent  flow  at  the  elevation  of  the  vent  as 
it  enters  the  top  space  from  the  bottom  space  and/or  as  it  enters  the  bottom 
space  from  the  top  space.  Rates  of  m.'-ss,  enthalpy,  and  products  of  combustion 
extracted  by  the  vent  flows  from  upp^  and  lower  layers  of  inside  room 

environments  and  from  outside  ambien  paces  are  determined  explicitly.  The 

algorithm/subroutine  is  called  VENTCF.  The  computer  subroutine  is  written  in 
FORTRAN  77.  The  subroutine  is  completely  modular,  and  it  is  suitable  for 
general  use  in  two -layer,  multi -room,  zone -type  fire  model  computer  codes.  It 
has  been  tested  over  a wide  range  of  input  variables  and  these  tests  are 
described. 


Keywords:  building  fires;  compartment  fires;  computer  models;  fire  models; 

mathematical  models ; vents ; zone  models 


1.  Introduction 


This  work  presents  an  algorithm  and  associated  computer  subroutine,  called 
VENTCF.  The  purpose  of  VENTCF  is  to  calculate  the  instantaneous  effects  on 
two- layer  fire  environments  of  the  quasi-steady  flow  through  a horizontal  vent 
connecting  two  spaces  involved  in  a compartment  fire.  The  subroutine  is 
designed  to  be  modular  and  easy  to  integrate  into  any  two -layer  zone -type 
compartment  fire  model. 

Flow  through  horizontal  vents  can  be  of  particular  importance  in  fires  on 
ships.  Because  of  the  generic  design  of  ships,  compartments  of  fire  origin 
can  be  connected  by  horizontal  vents  to  adjacent  spaces  on  a deck  below  and/or 
above.  As  examples  of  this  consider  the  case  of  a fire  in  an  otherwise 
poorly- ventilated  hold  of  a ship  during  the  staging  of  fire-fighting 
activities  initiated  through  a horizontal  vent  from  the  deck  above,  or  the 
case  of  a fire  in  a compartment  on  an  upper  deck  which  is  vented  through  an 
overhead  horizontal  penetration  to  the  outside  environment  on  the  top  deck. 
Knowledge  of  the  flow  exchanges  through  such  vents  are  important  since  the 
oxygen  convected  along  with  these  flows  can  be  critical  in  determining  the 
growth  (or  decay)  of  the  fire.  Also,  deck-to-deck  flows  through  such 
horizontal  vents  can  play  a major  role  in  determining  the  smoke  spread  and  the 
resulting  development  of  potentially  hazardous  conditions  thoughout  the  rest 
of  the  ship.  A major  motivation  for  this  work  is  to  provide  a computational 
tool  that  can  be  used  to  predict  these  phenomena. 

Consider  the  flow  through  a horizontal  ceiling/floor  vent  connecting  two 
spaces,  one  on  top  of  the  other,  involved  in  a fire  - generated  environment. 

The  two  spaces  can  be  two  inside  rooms  of  a multi-room  facility,  as  depicted 
in  Figure  la  of  the  APPENDIX,  or  they  can  be  comprised  of  one  inside  room  of  a 
facility  and  one  outside  space,  either  above  or  below  the  room,  which  is  used 
to  simulate  conditions,  local  to  the  vent,  of  the  outside  ambient.  The  latter 
configurations  are  depicted  in  Figures  lb  and  Ic  of  the  APPENDIX, 
respectively . 

The  flow  through  the  vent  is  determined  by  a model  which  takes  into  account 
traditional  considerations  of  flow  due  to  a cross -vent  pressure  difference 
(see,  e.g.,  reference  [1])  and,  when  appropriate,  flows  induced  when  existing 
fluid  configurations  are  unstable,  i.e.,  when  a relatively  cool,  dense  gas  in 
the  upper  space  overlays  an  elevated  temperature  and  relatively  less  dense  gas 
in  the  lower  space.  Whereas  a non-zero  cross-vent  pressure  difference  always 
tends  to  drive  a unidirectional  flow  through  the  vent,  from  the  higher-  to  the 
lower-pressure  side,  flow  induced  by  an  unstable  fluid  density  configuration 
tends  to  lead  to  an  exchange  - type  of  flow,  with  fluid  in  the  lower  space 
rising  into  the  upper  space  and  flow  from  the  upper  space  dropping  into  the 
lower  space.  In  the  present  algorithm,  the  actual  calculation  of  the  flow 
between  the  two  spaces  is  based  on  the  previously  developed  theory  presented 
in  reference  [ 2 ] . 

Determined  explicitly  by  the  present  algorithm/subroutine  are  the  net 
instantaneous  rates  of  addition  of  mass,  enthalpy,  and  products  of  combustion 
of  interest  and  the  properties  of  the  vent  flows  to  each  of  the  spaces  at  the 
elevation  of  the  vent.  A determination  of  where  such  flows  go  once  they  enter 
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the  recieving  spaces  would  be  obtained  with  the  use  of  additional  algorithms 
and  associated  subroutines.  These  additional  algorithms,  which  would  use  the 
output  of  the  present  algorithm/subroutine,  would  be  based  on  considerations 
beyond  the  scope  of  the  present  work. 

A flow  through  the  horizontal  vent  which  enters  one  of  the  spaces  joined  by 
the  vent  is  extracted  from  the  other  space.  Depending  on  the  configuration  of 
the  two  spaces,  the  direction  of  the  flow,  and,  in  the  case  of  inside  rooms, 
the  elevation  of  the  two-layer  interface  (i.e.,  at  the  floor,  ceiling,  or  in- 
between)  , the  present  algorithm  determines  explicitly  the  rates  of  ey'^action 
of  mass,  enthalpy,  and  products  of  combustion  from  the  upper  and  lowe  layers 
of  the  one  or  two  inside  rooms  and/or  from  an  outside  ambient  space  that  are 
joined  by  the  vent  under  consideration.  Along  with  other  components  of  flow 
to  or  from  the  layers  of  the  inside  rooms,  determined  with  the  use  of  other 
algorithms , these  rates  would  be  used  to  continue  in  time  the  solu:  jn  to  the 
equations  of  the  overall  fire  model.  These  are  the  equations  used  ^o  simulate 
mathematically  the  facility's  overall  dynamic  fire  environment. 


2.  The  Algorithm  and  Associated  Subroutine  VENTCF  for  Calculating  the  Effects 
of  Flow  Through  Horizontal  Ceiling/Floor  Vents 

The  algorithm/subroutine  VENTCF  is  presented  here  in  the  APPENDIX  as  a stand- 
alone document  that  can  be  inserted  as  a new  entry  into  the  catalog  of  modular 
algorithms  and  associated  computer  subroutines  of  reference  [3].  That 
document  is  a grouping  of  algorithms/subroutines  useful  for  simulating  the 
physical  phenomena  in  multi -room  zone -type  compartment  fire  model  computer 
codes . Pagination  of  the  present  APPENDIX  and  of  all  entries  of  the  reference 
[3]  catalog  is  according  to  name  of  the  particular  algorithm/subroutine  (in 
this  case,  VENTCF)  and  page  number  of  the  catalog  entry. 

The  catalog  of  reference  [3]  was  concieved  of  as  a growing  document  where  the 
entries  would  be  available  for  general  use  by  people  interested  in:  developing 
or  improving,  for  their  own  particular  needs,  a general  or  special-purpose 
multi -room  zone -type  compartment  fire  model;  or  predicting  isolated 
compartment  fire  phenomena,  for  whatever  reason.  The  present  VENTCF 
algorithm/subroutine  is  the  first  new  entry  to  the  original  reference- [ 3 ] 
document . 


3 . Format  of  the  VENTCF  Catalog  Entry 

The  development,  technology  transfer,  and  use  of  a reference -[ 3 ]- type  of 
catalog  'f  algorithms  and  cssociated  s ;broutines  is  enhanced  by  maintaining 
guidelines  for  a uniform  lormat  of  al  jrithm/subroutine  documentation.  In 
this  regard  a prototype  format  was  developed  and  used  in  all  reference- [ 3 ] 
algorithm/subroutine  catalog  entries.  The  format,  which  is  followed  here, 
includes  the  following  f ' ements : 
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TITLE 


TITLE 

Should  indicate  the  main  purpose  of  the 
algorithm/subroutine . 

DESCRIPTION 

General  description  of  the  algorithm. 

OUTPUT 

List  of  output  variables,  including  definitions 
and  units . 

INPUT 

List  of  input  variables,  including  definitions 
and  units . 

CALCULATIONS 

Concise  description  of  the  rules  for  obtaining 
the  output  variables  from  the  input  variables. 
This  would  include  or  refer  explicitly  to  all 
equations  required  in  the  calculation.  If  other 
algorithms/subroutines  are  required,  then  these 
should  be  readily  available  and  referenced. 

SUBROUTINES 

USED 

Listing  of  or  explicit  reference  to  each 
algorithm/subroutine  used  to  carry  out  the 
calculations . 

REFERENCES 

A list  of  references. 

SUBROUTINE 

VARIABLES 

A cross-reference  of  all  nomenclature  (including 
units)  introduced  in  the  above  sections  to  the 
nomenclature  used  in  the  FORTRAN  computer 
subroutine . 

PREPARED  BY 

Names  of  those  who  prepared  the  algorithm/sub- 
routine and  date  of  preparation. 

SUBROUTINE 

Listing  of  the  computer  subroutine.  This  would 
be  well -commented  and  would  include  a summary  of 
the  purpose  of  the  subroutine  and  definitions 
(including  units)  of  its  input  and  output 
variables. 

The  algorithm  VENTCF  and  a listing  of  its  associated  computer  subroutine  is 
presented  in  the  APPENDIX.  The  subroutine  is  coded  in  FORTRAN  77. 

4.  Testing  of  the  Subroutine  VENTCF 

The  subroutine  VENTCF  has  not  yet  been  incorporated  into  any  full  compartment 
fire  model  computer  code.  However,  extensive  parametric  testing  of  the 
subroutine  has  been  carried  out.  This  testing  will  be  described  here. 
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A wide  range  of  environment  scenarios  for  each  of  the  three  basic 
configurations  of  Figure  1 of  the  Appendix  were  considered.  For  each 
configuration,  parameters  were  varied  in  a manner  as  to  simulate  all  possible 
combinations  of  the  following: 

a.  for  the  two  layers  of  an  inside  space:  the  usual  case  with  two 
non- zero  - thickness  layers  (i.e.,  layer  interface  between  the 
ceiling  and  the  floor) , or  one  non-zero- thickness  layer  and  one 
zero  - thickness  layer  (i.e.,  layer  interface  at  the  ceiling  or 
floor) ; 

b.  stable  or  unstable  cross -vent  density  configuration; 

c.  the  reference  elevation  for  an  outside  space  is  above,  at,  or 
below  the  vent  elevation; 

d.  1.0  atmosphere  reference  pressure,  p^^p  space  (i.e., 

Pdatum  ^ 101325.  pa  with  5pj^pp  = Pr£p  i ■ Pdatum  with 

reference  pressure  in  the  bottom  space,  p^pp  2>  varying  from  0.01 
atmospheres  to  2.0  atmospheres  (i.e.,  < 5pj^pp  2 — 

Pdatum  ^ 

e.  1.0  atmosphere  reference  pressure,  P^pp  £>  bottom  space 

(i.e.,  Pdatum  ^ 101325.  pa  with  5pj^pp  2 - Pref,2  " Pdatum  ^ 
with  reference  pressure  in  the  top  space,  p^pp  , varying  from 
0.01  atmospheres  to  2.0  atmospheres  (i.e.,  "O-^^Pq^^um  ^ ^Pref  1 - 
Pdatum  ) 

Five  distinct  elevations  (relative  to  a datum  elevation),  yj^  , N = 1 to  5 , and 
three  distinct  densities,  , N = 1 to  3 , were  required  to  construct  the 
environment  scenarios.  These  were  chosen  to  be: 


y^  = 0.0m;  y2  = 1.5m;  y3  = 3.0m;  y^  = 4.5m;  y^  = 6 . Om 
p^  = l.OOkg/m^;  = 0.50kg/m^;  p^  = 0.25kg/m^ 


(1) 


When  the  environment  of  an  inside  room  involved  two  layers,  the  layer 
densities  were  specified  to  be  in  a stable  configuration  (i.e.,  a low  density 
layer  above  a high  density  layer) . 

Table  1 identifies  all  configurations  and  environment  specifications  that  were 
used  in  the  parametric  testing  of  the  subroutine.  This  involved  162  basic 
"Cases."  Each  Case  involved  a separate  parametric  study  of  over  two  thousand 
calls  to  the  subroutine.  These  covered  the  above  pressure  range  of  either 
condition  d.  or  condition  e.  The  calls  involved  incremental  changes  of  the 
relative  reference  pressures,  <5pj^pp  for  condition  e and  5pj^pp  2 
condition  d,  which  were  small  enough  to  reveal  details  of  the  exchange- flow 
phenomena,  these  being  sensitive  to  fractional -pascal- level  variations  in 
cross-vent  pressure  differences. 
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For  each  Case,  Table  1 identifies  values  of  the  parameters  which  define  the 
configuration  of  the  two  spaces,  as  illustrated  in  Figure  1 of  the  Appendix, 
and  the  states  of  the  environments  in  the  spaces.  These  are  the  parameters 
used  to  call  the  subroutine.  (The  reader  is  referred  to  the  NOMENCLATURE 
Section  for  explaination  of  terms  used  in  the  table.)  Also  explicitly- 
identified  in  the  table  are  those  Cases  which  involve  unstable  cross -vent 
density  configurations.  When  cross-vent  pressure  differences  are  small  enough 
in  magnitude,  such  Cases  lead  to  the  exchange  - flow  phenomena  which  are  of 
particular  interest. 

Presented  in  Figures  1 and  2 are  plots  of  the  calculated  mass  flow  rates  to 
upper  room  1,  and  to  lower  room  2,  as  functions  of  the 

computed  cross -vent  pressure  difference,  Ap , for  Cases  28,  29,  and  30.  As  in 
Eq.  (2)  of  the  Appendix,  Ap , an  output  variable  of  the  subroutine,  is  defined 
as 


Ap  = P2  - Pi 


(2) 


where  p^  and  pressures  in  the  top  and  bottom  spaces  at  the 

elevation  of  the  vent  (see  Figure  2 of  the  Appendix). 

As  indicated  Table  1,  Cases  28,  29,  and  30  involve  two  inside  rooms  where,  as 
depicted  in  Figure  la  of  the  Appendix,  both  of  these  have  an  upper  and  a lower 
layer.  Case  28  involves  a neutrally- stable  cross -vent  density  configuration 
since  the  density  immediately  above  the  vent  (in  the  lower  layer  of  room  1), 

Pl  i ^ Pi  ^ 0.50kg/m^,  is  identical  to  the  density  immediately  below  the  vent 
(in  the  upper  layer  of  room  2).  Cases  29  and  30  involve  an  unstable  cross- 
vent density  configuration  since  the  density  immediately  above  the  vent,  ^ 

= Pi  = 1 ■ Okg/m^ , is  greater  than  the  density  immediately  below  the  vent,  p^  ^ 

= P2  = 0 • 50kg/m^ . 

The  scale  of  the  Ap  abscissa  of  Figure  1 is  relatively  coarse,  and  the  plots 
reveal  the  calculated  mass  flow  rates  through  the  vent  as  a result  of  large 
cross -vent  pressure  differences,  i.e.,  when  effects  of  compressibility  become 
significant  as  could  be  the  case,  for  example,  for  fires  in  nearly- 
hermetically- sealed  facilities.  In  the  figure,  the  sharp  breaks  in  the 
^VENT  1 ^VENT  2 P^o^s  for  Cases  28  and  30  and  for  Case  29,  respectively, 

correspond  to  cross -vent  pressure  differences,  Ap , which  separate  choked  from 
unchoked  flow  through  the  vent.  Define  Phiqh  and  Plow  absolute 

pressures  on  the  high-  and  low-pressure  side  of  the  vent,  respectively.  Then, 
consistent  with  Eqs . (7) -(15)  of  the  Appendix,  such  choking  occurs  for  air 
(i.e.,  7 = 1.40)  when  [4,5] 

Plow/Phigh  ^ [2/(7  + I)]’''’  - - 0.528  (3) 


As  indicated  in  Figure  1,  for  Case  29  the  parameters  specified  in  Table  1 lead 
to  a prediction  of  choked  vent  flow  from  room  1 to  room  2 approximately  when 
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Ap  < - 0 . 472pq^^  uj,  -=  -0 . 478  ( 10^  )pa . Also,  for  Cases  28  and  30,  choked  flow 
from  room  2 to  rc  1 is  predicted  when  Ap  > 0 = 0 . 478  ( 10^  )pa . 

The  scale  of  the  oscissa  of  Figure  2 is  relatively  fine,  and  the  plots 

reveal  the  calcula_  . mass  flow  rates  through  the  vent  when  the  cross -vent 
pressure  differences  are  at  or  close  to  zero.  Since  the  Case-28  density 
configuration  is  neutrally- stable , i.e.,  only  uni-directional  or  zero  flow  is 
possible,  the  mass  flow  rates  through  the  vent  corresponding  to  Ap  = 0 are 
seen  to  be  identically  zero.  However,  for  the  unstable  cross -vent  density 
configurations  of  Cases  29  and  30,  the  vent  flow  algorithm  is  seen  t d to 

the  exchange  - flow  phenomenon.  In  particular,  there  is  a non- zero  ve  ow 

exchange  between  the  top  and  bottom  spaces  whenever  | Ap  | < App^ooD  ’ v--- - - 
AppLooD  defined  in  Eq . (19)  of  the  Appendix. 


ACKNOWLEGEMENTS 

The  author  acknowleges  gratefully  the  Naval  Research  Laboratory  which 
sponsored  this  work  under  the  technical  management  of  Dr.  Patricia  A.  Tatem 
and  with  the  support  of  Ms.  Jean  L.  Bailey. 


REFERENCES 


[1]  Enunons , H. , SFPE  Handbook  of  Fire  Protection  Engineering.  Sect. 

1/Chapter  8,  SFPE,  Boston,  pp . 130-138,  1988. 

[2]  Cooper,  L.Y. , Calculation  of  the  Flow  Through  a Horizontal  Ceiling/Floor 
Vent,  NISTIR  89-4052,  National  Institute  of  Standards  and  Technology, 
Gaithersburg  MD , March  1989. 

[3]  Cooper,  L.Y.  and  Forney,  G.P.,  Editors,  The  Consolidated  Compartment 
Fire  Model  (CCFM)  Computer  Code  Application  CCFM. VENTS  - Part  III: 
Catalog  of  Algorithms  and  Subroutines,  NISTIR  90-4342,  National 
Institute  of  Standards  and  Technology,  Gaithersburg  MD,  1990. 

[4]  Cooler,  L.Y. , Calculating  Flows  Through  Vertical  Vents  in  Zone  Fire 
Moce.;,  Combustion  Science  and  Technology,  Vol.  63,  Nos.  1-3,  pp . 43-50, 
1989. 

[5]  Shapiro,  A.H. , The  Dynamics  and  Thermodynamics  of  Compressible  Fluid 
Flow.  Rowland  Press,  1953. 


7 


NOMENCLATURE 


M 

'•^VENT  , 


Pd ATUM 


Ph  I GH 


Plow 


Pref  , I 


Pi  [Pz 


Yc-  Yl 


yCE  IL  . 


Flayer 


j ; I = 1 to  2 

Mass  flow  rate  of  the  vent  flow  component  entering  space  I 


Datum  absolute  pressure 


Maximum  of  (p^ , P2 ) 


Minimum  of  (p^ , P2 ) 

; I = 1 or  2 

Absolute  pressure  in  space  I at  the  reference  elevation  y^gp  j 


Absolute  hydrostatic  pressure  in  the  top  [bottom]  space  at  the  elevation 
of  the  vent 


, for  space  I 

ycEiL.i’  ytAYER,!-  ^^^yREF,!’  respectively. 

I [yREF , I ^ ’ I = 1 or  2 

If  space  I is  an  inside  room:  elevation  of  the  ceiling  [floor]  of  room 
I above  the  datum  elevation;  if  space  I is  an  "outside  room":  ycEiL  i 
and  y^pp  j are  both  identical  and  equal  to  the  reference  elevation  of 
space  I above  the  datum  elevation,  i.e.,  y^piL  i ^ yREF  i 


j ; I = 1 or  2 

If  space  I is  an  inside  room:  elevation  of  the  upper/lower  layer 
interface  in  room  I above  the  datum  elevation;  if  space  I is  an  "outside 
room":  ypAYER,!  ^ yREF,I 
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YvE  N T 


Elevation  of  the  vent  above  the  datum  elevation 


Ap 


Pz  - Pi 


^Pf  lood 

Minimum  v e of  | Ap | , in  cases  of  unstable  cross-vent  density 
conf  igurat i.ons , leading  to  unidirectional  vent  flow;  Eq . (19)  of 

Appendix 


^Pref , I ’ ^ 1 or  2 

Pressure  at  the  reference  elevation,  j , in  space  I above  the  datum 

absolute  pressure,  Pd^tum  ’ space  I is  an  inside  room,  then  5pj^j.p  j 
and  y^Ep  j correspond  to  the  pressure  and  elevation,  respectively, 
within  the  room  and  at  the  floor 


Pl  , I t Pu  , I 1 ^ ^ ^ 

If  space  I is  an  inside  room:  density  of  the  lower  [upper]  layer  in 
room  I if  the  volume  of  the  lower  [upper]  layer  is  non-zero  (if  the 
lower  [upper]  layer  volume  is  zero  then  the  Pl  i value  is  not 

used  in  the  calculation);  if  space  I is  an  "outside  room":  j is  the 

uniform  density  there  and  Py  i = Pl  i 
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TABLE  1: 


Conf igurations/Environments  Used  to  Test  the  VENTCF  Subroutine 


= l.OOkg/m^  > P2  0.50kg/m^  > Pg  = 0.25kg/m^ 
y,  =O.Om<y,  =1.5m<y-  =3.0m<y,  =4.5m<y^  =6.0m 
pL.uh  - 1.0i!25(10=)pa 


CASES  1-54:  SPACE  1 (TOP  SPACE)  IS  INSIDE  ROOM,  SPACE  2 (BOTTOM  SPACE) 

IS  INSIDE  ROOM  - APPENDIX:  FIGURE  la 


CASE  SPACE  SPACE 
NO.  1 2 


1 

2 

3 

4 

5 

6 


7 

8 

9 

10 
11 
12 


13 

14 

15 

16 

17 

18 


Lower  Lower 
Layer  Layer 


yR=y3  yR=yi 
yL=y5  yL=y3 
yc=y5  yc=y3 


Pl.  1=P2<Pl.2=Pi 

^Pref , • 

» "Pdatum*^^Pref  , 2— Pdatum 

n 

^Pref  , 2""^  • 

» “Pdatum'^^Pref  , 1— Pdatum 

Pl,  i=Pi=Pl.2 

^Pref , • 

» "Pdatum^^Pref , 2-Pdatum 

ft 

^Pref  , 2""^  • 

» "Pdatum^^Pref , 1— Pdatum 

Pl.i=Pi>Pl.2=P2 

^Pref  , • 

> “Pdatum^^Pref , 2— Pdatum 

n 

^PrEF , 2^^ • 

» "Pdatum'^^Pref  , 1-Pdatum 

* 

* 


Two 

Layers 


yR=yi 

yL=y2 

yc=y3 


Pl.  ,2  Pz^P'L.Z  Pi 

ft 


Pl  , l~P2’^Pu  , 2^Pl  , 2 Pi 

n 


P2^Pl.2  Pi 
i 


^PREF  , 1""®  • ’ “PdATUM'^^PrEF  , 2— PdATUM 

^Pref  , 2^^'  ’ “Pdatum^^Pref  , 1— Pdatum 
^pREF  , * » "Pd  ATDM'^^PrEF  , 2— Pd  ATOM 

^Pref  , 2^^'  » "Pdatom'^^Pref  , 1-PdATUM 
^Pref  , i"'®  • » "Pdatum'^^Pref  , 2— Pdatum 
^Pref  , 2^*  » "Pdatum'^^Pref  , 1— Pdatum 


Upper 

Layer 

yR=yi 

yL=yi 

yc=y3 


Pl,  l=P2<Pu,2=Pl 

^Pref , • 

» "Pdatum'^^Pref  , 2— Pdatum 

t1 

^PrEF ,2^^' 

» "Pdatum'^^Pref  , 1-Pdatum 

Pl.  l=Pl=Pu.2 

^Pref  . 1""®  * 

» "Pdatum'^^Pref  , 2-Pdatum 

n 

^Pref  . 2""®  • 

» “Pdatum“^^Pref  , 1-Pdatum 

Pl  , l^Pl^Pu  , 2^p2 

^Pref  . 1""®  • 

» "Pdatum'^^Pref  , 2-Pdatum 

n 

^Pref  , 2""®  • 

» "Pdatum'^^Pref  , 1— Pdatum 

* Unstable  cross -vent  density  configuration 
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TABLE  1:  Conf igurations/Enviroiunents  Used  to  Test  the  VENTCF  Subroutine 

(Cont ' d) 

= l.OOkg/m^  ^ Pz  ^ 0.50kg/m^  > Pg  = 0.25kg/m^ 
y,  = 0 . Om  < y,  = 1 . 5m  < y-,  = 3 . Om  < y*  = 4 . 5m  < y-  = 6 . Om 
= 1.0ll25(10=)pa 


CASES  1-54:  SPACE  1 (TOP  SPACE)  IS  INSIDE  ROOM,  SPACE  2 (BOTTOM  SPACE) 

(CONT'D)  IS  INSIDE  ROOM  - APPENDIX:  FIGURE  la 


CASE  SPACE  SPACE 
NO.  I 2 


Two  Lower 
Layer  Layer 

yR=y3  yR=yi 
yL=y4  yL=y3 
yc=y5  yc=y3 


19 

Pl , 1 Pz^Pl , 2 Pl 

20 

tl 

21 

Pl  , i^Pi^Pl  , 2 

22 

n 

23  * 

Pl  , i^Pi^Pl  . z^Pz 

24  * 

4. 

^Pref ,1  ^ » Pd  ATUM  “^^Pref  , 2— Pd  ATUM 
^PrEF  , 2^  • ’ “PdATUM'^^PrEF  , 1— Pd  ATUM 
^pREF  , • ’ “PdATUM"^^PrEF  , 2” PdATUM 

^Pref  , 2"^^  • » "Pdatdm'^^Pref  , i— Pdatum 
^Pref  . • ’ “Pdatum^^Pref  , 2 -Pd  at  dm 

^Pref  , 2^®  • » 'Pdatum^^Pref  , 1— Pdatum 


25 

26 

27 

28 

29 

30 


31 

32 

33 

34 

35 

36 


* 

* 


* 

* 4. 


Two 

Layer 

yR=yi 

yL=y2 

yc=y3 


4. 


P'L  . i^P2^P\} . z^Pz^Pl  , z^Pi  ^Pref  , 1""^  • » "Pdatum'^^Pref  , 2-Pdatdm 

^pREF  . 2""^  • » "Pdatum'^^Pref  , 1-Pdatdm 
Pl  , l^Pz^Pu  . Z^Pl.  . Z^Pl  ^pREF  , • • “PdATDm'^^PrEF  , 2“ PdATDM 

^Pref  , 2"^' » "Pdatdm“^^Pref  , 1-Pdatdm 
Pl  , i""Pi^Pu  , z^Pz^Pl  . z^Pi  ^Pref  , » "Pdatdm^^Pref  , 2-Pdatdm 

^pREF  , 2""^  • » "Pdatdm'^^Pref  , 1-Pdatum 


Upper 

Layer 

yR=yi 

yL=yi 

yc=y3 


Pl  , 1 Pz^Pu , 2 Pl 

^Pref , 1^^  * 

» "Pdatdm'^^Pref  , 2— Pdatum 

n 

^Pref , 2^^ • 

» "Pdatum“^^Pref  , 1— Pdatum 

Pl , i^Pi^Pu , 2 

^Pref , 1^^ • 

» "Pdatdm'^^Pref  . 2-Pdatdm 

tl 

^Pref  , 2""^  • 

’ ■Pdatum'^^Pref  , 1— Pdatum 

Pl  . i^Pi^Pu  . z^Pz 

^Pref.i=0* 

» "Pdatum^^Pref , 2-Pdatdm 

n 

^Pref  , 2""^  ‘ 

» "Pdatdm'^^Pref  , 1— Pdatum 
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TABLE  1:  Conf igurations/Environments  Used  to  Test  the  VENTCF  Subroutine 

(Cont ' d) 


= l.OOkg/m^  ^ Pz  0.50kg/m^  > Pg  = 0.25kg/m^ 

y,  = 0.0m  < y,  = 1.5m  < y-,  = 3.0m  < y.  = 4.5m  < y-  = 6.0m 
- l.of325(10=)pa 


CASES  1-54:  SPACE  1 (TOP  SPACE)  IS  INSIDE  ROOM,  SPACE  2 (BOTTOM  SPACE) 

(CONT'D)  IS  INSIDE  ROOM  - APPENDIX:  FIGURE  la 


CASE  SPACE  SPACE 
NO.  1 2 


37 

38 

39 

40 

41 

42 


43 

44 

45 

46 

47 

48 


Upper  Lower 
Layer  Layer 


yR=y3  yR=yi 
yL=y3  yL=y3 
yc=y5  yc=y3 


Pu  , 1 Pz^Pl  . 2 Pi 

^Pref  , i""®  • 

» "Pdatum'^^Pref  , 2“PdATUM 

n 

^PREF  , 2^®  • 

’ “Pdatum^^Pref , 1“Pdatum 

P\J  . i'^Pi'^Pl  . 2 

^Pref  . 1=0- 

» "Pdatum^^Pref , 2—Pdatum 

tt 

^PrEF , 2^® • 

» "Pdatum'^^Pref  , 1— Pdatum 

P\j  , I^Pi^Pl  , Z^Pz 

^Pref , • 

» “Pdatum'^^Pref  , 2—Pdatum 

n 

^PrEF , 2^^ • 

» "Pdatum'^^Pref  , 1— Pdatum 

* 

* 


Two 

Layer 


yR=yi 

yL=y2 

yc=y3 


i 


, i^Pz^Pv  , z^Pz'^Pl  . z^Pi  ^Pref  , 1°"^  * ’ "Pdatum'^^Pref  , 2— Pdatum 

^Pref  . 2^®*  » "Pdatum'^^Pref  . 1-Pdatdm 
P\j , , 2'^Pl  , z^Pi  ^Pref  , * ’ "Pdatum'^^Pref  , 2— Pdatum 

^pREF  , 2“®  • » "PdATUM'^^PrEF  , I-PdATUM 

P\j  , l^Pl^Pv  . Z^Pz^Pl  . Z^Pl  ^PREF  . 1^^'  » "PdATUM'^^PrEF  , 2-PdATUM 

^Pref  , 2""®  • ’ "Pdatum'^^Pref  , 1-Pdatdm 


49 

50 

51 

52 

53 

54 


Upper 

Layer 


yR=yi 

yL=yi 

yc=y3 


P\J  . 1 Pz^^Pu  . 2 Pi 

ft 


, 1 Pi  P\J  , z 

tt 

* Pu  , 1 ^Pl^P\J  , 2 ^Pz 

* i i " 


^Pref  , • » "Pdatum'^^Pref  , 2— Pdatum 
^Pref  . 2^®  • ’ "Pdatum'^^Pref  , 1 “Pdatum 
^Pref  , * ' "Pdatum'^^Pref  , 2—Pdatum 

^Pref  , 2"^^  * ’ "Pdatum'^^Pref  , 1— Pdatum 
^Pref  , • ’ “Pdatum'^^Pref  , 2—Pdatum 

^Pref  , 2^^ • » "Pdatum'^^Pref  , 1“ Pdatum 
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TABLE  1:  Conf igurations/Environments  Used  to  Test  the  VENTCF  Subroutine 

(Cont ' d) 


= l.OOkg/m^  Pz  ^ 0.50kg/m^  > Pg  = 0.25kg/m^ 

Vt  = 0 . Om  < y-  = 1.5m  < = 3 . Om  < y . = 4.5m  < y-  = 6,0m 

Pdatuh  = 1.01325(105)pa 


CASES  55-108:  SPACE  1 (TOP  SPACE)  IS  OUTSIDE  SPACE,  SPACE  2 (BOTTOM  SPACE) 

IS  INSIDE  ROOM  - APPENDIX:  FIGURE  lb 


CASE  SPACE  SPACE 
NO.  1 2 


All  Lower 
One  Layer 
Layer 


yR=y4  yR=yi 
yL=y4  yL=y3 
yc=y4  yc=y3 


55 

Pl  . 1 Pz^P\.  , 2 Pi 

56 

n 

57 

Pl  , I^Pi^Pl  , 2 

58 

n 

59  * 

Pl  , i^Pi^Pl  . z^^Pz 

60  * 

n 

Two 

Layers 

yR= 

=yi 

yL= 

=y2 

yc= 

=y3 

61 

Pl  , i~P3“^Pu  , z^Pz'^Ph  , 2^P 

62 

n 

63 

Pl  . i""P2""Pd  , 2'^Pl  , 2^Pl 

64 

n 

65  * 

P-L  . l~Pl^Pu  . Z~P2^Pl  , 2^P 

66  * 

4. 

Upper 

Layer 

yR^yi 

yL= 

=yi 

yc^ 

-y3 

67 

Pl  , i'"P2"^Pu  , 2""Pi 

68 

n 

69 

Pl  , i^Pi""Pu  , 2 

70 

n 

71  * 

Pl  , I'^Pi^Pu  , 2^P2 

12  * i 4. 


^Pref  , i""®  • » "Pd atum'^^Pref  , 2— Pdatum 
^Pref  , 2^®  * » "Pdatdm'^^Pref  , 1— Pdatum 
^Pref  , • » "Pd  ATUM  ^^Pref  , 2— Pda  TDM 
^Pref  , 2""^  • » "Pdatum'^^Pref  , 1— Pdatum 
^Pref  , • » "Pdatdm^^Pref  , 2— Pdatum 
^Pref  , 2"^  • ’ "Pdatum"^^Pref  , 1— Pdatum 


^Pref  , • » "Pdatum'^^Pref  , 2— Pd atdm 
^Pref  , 2"^  • ’ "Pdatum'^^Pref  , 1— Pdatum 
^Pref  , * » "Pdatum"^^Pref  , 2— Pdatum 
^Pref  , 2^  • ’ "Pdatdm'^^Pref  , 1-Pdatdm 
^Pref  , • » "Pdatum'^^Pref  , 2-Pdatdm 

^Pref  , 2^  • ’ "Pdatdm'^^Pref  , 1-Pdatdm 


^Pref  . 1^-  » "Pdatum'^^Pref  , 2-PdATDM 
^Pref  , 2^^  • » "Pdatdm^^Pref  , 1— Pdatum 
^Pref  , • » "Pdatdm'^^Pref  , 2— Pdatum 

^Pref  , 2^  • ’ "Pdatdm'^^Pref  , 1-Pdatdm 
^Pref  , * ’ "Pdatum'^^Pref  , 2— Pdatum 

^Pref  , 2^^  ■ ’ "Pdatdm'^^Pref  , 1— Pdatum 
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TABLE  1:  Configurations/Enviroraoents  Used  to  Test  the  VENTCF  Subroutine 

(Cont^d) 

= l.OOkg/m^  > Pi  0.50kg/m^  > p^  = 0.25kg/m^ 

y,  = 0.0m  < ==  1.5m  < y-  = 3.0m  < y,  = 4.5m  < y^  = 6.0m 

= 1.0a25(10=)pa 


CASES  55-108:  SPACE  1 (TOP  SPACE)  IS  OUTSIDE  SPACE,  SPACE  2 (BOTTOM  SPACE) 

(CONT'D)  IS  INSIDE  ROOM  - APPENDIX:  FIGURE  lb 


CASE  SPACE  SPACE 
NO.  1 2 


All  Lower 
One  Layer 
Layer 

yR=y3  yR=yi 


yL= 

=y3  yL= 

=y3 

yc= 

=y3  yc= 

=y3 

73 

Pl.  i=P2<Pl.2=Pi 

74 

ft 

75 

Pl  . . 2 

76 

n 

77  * 

^l.i=Pi>Pl.2=P2 

00 

4. 

Two 

Layer 

yR= 

=yi 

yL= 

-Vz 

yc= 

=y3 

79 

Pl  , l^P3^Pu  , 2~Pz^Pl  , Z^Pl 

80 

n 

81 

Pl  , I^Pz^Pm  , 2'^Pl  , 2^Pl 

82 

n 

83  * 

Pl  . i""Pi^Pu  . z^Pz^Pl  . z^Pi 

84  * 

i 

Upper 

Layer 

yR= 

“yi 

yL= 

“yi 

yc-y3 

85 

Pl  , l~P2'^Pu  . 2 Pi 

86 

n 

87 

Pl  , i^'Pi^Pu  , 2 

88 

n 

89  * 

Pl  , i""Pi^Pu  , z""Pz 

90  * i 4. 


^PREF ,1  ^ ’ Pd  ATUM  “^^Pref  , 2— Pd  ATUM 

^Pref  , 2""^  • » “Pd  atum^^Pref  1— Pdatum 
^Pref , • » “Pdatum^^Pref , 2— Pdatum 
^Pref , 2^® • » “Pdatum^^Pref , 1— Pdatum 
^Pref , • ’ "Pdatum^^Pref , 2— Pdatum 
^Pref  , 2^®  * » “Pdatum'^^Pref  , 1— Pdatum 


^Pref  , • » “Pdatum'^^Pref  , 2— Pdatum 
^Pref  , 2^®  • » “Pdatum'^^Pref  , 1— Pdatum 
^Pref  , • ’ “Pdatum^^Pref  , 2— Pdatum 
^Pref  , 2""®-  » “Pdatum'^^Pref  , 1— Pdatum 
^Pref  , i""®  • » “Pdatum'^^Pref  . 2— Pdatum 
^Pref  , 2^®*  » “Pdatum'^^Pref  , 1— Pdatum 


^Pref  , * » “Pdatum'^^Pref  , 2— Pdatum 

^Pref  , 2"“®  • » “Pdatum'^^Pref  , 1— Pdatum 
^Pref  . ’ “Pdatum*^^Pref  , 2-Pdatum 

^Pref  , 2""®*  » "Pdatum'^^Pref  , 1-Pdatum 
^Pref  , • » “Pdatum'^^Pref  , 2— Pdatum 

^Pref  , 2^®  * » “Pdatum'^^Pref  , 1— Pdatum 
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TABLE  1:  Conf igurations/Environments  Used  to  Test  the  VENTCF  Subroutine 

(Cont ' d) 


= l.OOkg/m^  ^ Pt.  ""  0.50kg/m^  > p^  = 0.25kg/m^ 
y,  = 0 . Om  < y,  = 1 . 5m  < y,  = 3 . Om  < y , = 4 , 5m  < y..  = 6 . Om 
- 1.0l!25(10=)pa 


CASES  55  08:  SPACE  1 (TOP  SPACE)  IS  OUTSIDE  SPACE,  SPACE  2 (BOTTOM  SPACE) 

(CONT'D)  IS  INSIDE  ROOM  - APPENDIX:  FIGURE  lb 


CASE  SPACE  SPACE 
NO.  1 2 


91 

92 

93 

94 

95 

96 


97 

98 

99 

100 
101 
102 


103 

104 

105 

106 

107 

108 


All  Lower 
One  Layer 
Layer 


yR==y2  yR=yi 
yL=y2  yL=y3 
yc=y2  yc=y3 


Pu  .i==P2<Pl,2=Pi 

tt 


. 1 Pi  Pl .2 

n 

P\J  , I^Pi^Pl  , 2^P2 

n 


* 

* i 


Two 

Layer 

yR=yi 

yL=y2 

yc=y3 


, l“^3‘^^U  , 2 P2'^Pl,2  Pi 

n 

P\J  . l^P2^PlS  . 2"^Pl  . 2^Pl 

tf 

Pm  , I^Pi^Pm  , 2^P2^Pl  . 2^Pl 


Upper 

Layer 

yR=yi 

yL=yi 

yc='y3 


P\J  . 1 P2^P\J  , 2 Pl 

n 

, l^Pl""P\J  , 2 

P\J  . l""Pl^Pu  . 2^P2 


^Pref  , • 

» "Pdat 

^PrEF , 2^^ • 

» "Pdat 

^Pref , • 

» "Pdat 

^PrEF ,2^^- 

' "Pdat 

^Pref , * 

’ "Pdat 

^PrEF . 2^^ • 

’ "Pdat 

um'^^Pref  , 2-Pdatum 
dm^^Pref  , 1— Pdatum 
um'^^Pref  , 2-PdATUM 
UM  <«Pref.i^Pda  TUM 
dm"^^Pref  , 2— Pdatum 
um'^^Pref  , 1— Pdatum 


^Pref  , • » “Pdatum"^^Pref  . 2— Pdat 

^Pref  , 2^®  • » "Pdatum'^^Pref  , 1— Pdat'u 
^Pref  , * ’ "Pdatum"^^Pref  , 2— Pdatu*-: 

^Pref  , 2^®  • » "Pd atum'^^Pref  , 1— Pdatum 
^Pref  , • » "Pdatum^^Pref  , 2— Pdatum 

^Pref  , 2""^  • ’ "Pdatum^^Pref  , 1— Pdatum 


^Pref  , • ’ “Pdatum^^Pref  , 2— Pdatum 

^Pref  ,2^^'  » "Pdatum'^^Pref  , 1— Pdatum 

^pREF  , * » "PdATUM'^^PrEF  , 2“ PdATUM 

^Pref  , 2^^  • ’ “Pdatum'^^Pref  , 1— Pdatum 
^Pref  , • » "Pdatum'^^Pref  , 2— Pdatum 
^Pref  , 2^^  • » "Pdatum“^^Pref  , 1— Pdatum 
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TABLE  1;  Configurations /Environments  Used  to  Test  the  VENTCF  Subroutine 

(Cont ' d) 

Pi  = l.OOkg/m^  > ^2  ""  0.50kg/m^  > pg  = 0.25kg/m^ 
y,  =O.Om<y,  =1.5m<y,  =3.0m<y.  =4.5m<y-  =6.0m 
pLtum  = 1.0i325(10=)pa 


CASES  109-162:  SPACE  1 (TOP  SPACE)  IS  INSIDE  ROOM,  SPACE  2 (BOTTOM  SPACE) 

IS  OUTSIDE  SPACE  - APPENDIX:  FIGURE  Ic 


CASE  SPACE  SPACE 
NO.  1 2 


109 

110 
111 
112 

113 

114 


Lower  All 
Layer  One 

Layer 


yR=y3  yR=y4 
yL=y5  yL=y4 
yc=y5  yc=y4 


P-L  .2  Pi 

tt 


* 

* 


, 1 Pi  P\.  ,2 

t1 

, l^Pl^Ph  . 2^P2 

4. 


^Pref  , * 

> “Pdat 

^Pref , 2^® • 

» "Pdat 

^Pref  , 1""®  • 

» “Pdat 

^Pref  , 2^^  • 

» “Pdat 

^Pref  , 1""®  • 

» “Pdat 

^Pref , 2^® • 

» “Pdat 

um'*^^Pref  . 2— Pdatum 
um'^^Pref  , 1— Pdatum 
UM  “^^Pref  , 2— Pdat  UM 
um'^^Pref  , 1— Pdatum 
um'^^Pref  . 2— Pdatum 
um'^^Pref  , 1-Pdatum 


115 

116 

117 

118 

119  * 

120  * 


yR=y3 

yL=y3 

yc=y3 


i 


Pl  , 1 Pz^Pm.Z  P2^P\.,2  Pi 

tt 

Pl  , l^P2""Pu  , 2^Pl.  , 2^Pl 

n 

Pl  . l""Pl^Pu  . Z^Pz^Pl  . Z^Pl 


^Pref  ,1^®-  » "Pdatum'^^Pref  , 2— Pdatum 
^pREF  . 2""^*  » "PdATUM'^^PrEF  , I" PdATUM 
^Pref  , i""®  • » "Pdatum'^^Pref  , 2— Pdatum 
^Pref  , 2""®-  » "Pdatum'^^Pref  , 1— Pdatum 
^Pref  , 1""®  • » "Pdatum'^^Pref  , 2-Pdatum 
^Pref  . 2^®  • » “Pdatum"^^Pref  . 1— Pdatum 


121 

122 

123 

124 

125  * 

126  * i 


yR=y2 

yL=y2 

yc=y2 


Pl  . i~Pz^P\j  . 2~Pi 

n 


Pl.i^Pi^Pm.z 

tt 

Pl  . i^Pi^P\5  . z^Pz 
i 


^Pref  , 1^®  • » "Pdatum'^^^Pref  , 2— Pdatum 
^Pref  . 2^  • » "Pdatum'^^Pref  , i-Pdatum 
^Pref  , 1^®*  » 'Pdatum'^^Pref  , 2-Pdatum 
^Pref  , 2^*  » "Pdatum^^Pref  . 1-Pdatum 
^Pref  , 1"^®  • » "Pdatum'^^Pref  , 2— Pdatum 
^Pref  , 2"^  • » "Pdatum'*^^Pref  , 1— Pdatum 
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TABLE  1:  Conf igurations/Environments  Used  to  Test  the  VENTCF  Subroutine 

(Cont ' d) 

= l.OOkg/m^  ^ Pz  ^ 0.50kg/m^  > pg  = 0.25kg/m^ 
y.  =0.0m<y,  =1.5m<yo  =3.0in<y.  = 4.5in<y.-  =6.0m 
= 1.0f325(10=)pa 


CASES  109-162:  SPACE  1 (TOP  SPACE)  IS  INSIDE  ROOM,  SPACE  2 (BOTTOM  SPACE) 

(CONT'D)  IS  OUTSIDE  SPACE  - APPENDIX:  FIGURE  Ic 


CASE  SPACE  SPACE 
NO.  1 1 


127 

128 

129 

130 

131 

132 


Two  All 
Layer  One 

Layer 

yR=y3  yR=y4 
yL=y4  yL=y4 
yc=y5  yc=y4 


P-L  , l~Pz^P'L  . 2 Pi 

n 


Pi  .1  Pi  Pi  .2 

n 


* 

* 


Pi  . 1 Pl^Pl  . 2 Pi 
i 


^PrEF  , l""®  • » "PdATUM'^^PrEF  , 2-PDATaM 
^Pref  , 2""^*  » "Pdatum'^^Pref  , 1— Pdatum 
^Pref  , i""®  • » "Pdatum'^^Pref  , 2— Pdatdm 
^Pref  , 2^^  • » "Pdatum'^^Pref  , 1— Pdatum 
^Pref  . • » "Pdatum'^^Pref  , 2— Pdatum 

^pREF  , 2^®  • » "Pdatum'^^Pref  , I” Pdatum 


133 

134 

135 

136 

137  * 

138  * 


yR=y3 

yL=y3 

yc^y3 


i 


Pi  , 1 Pz^Pv.Z  Pl'^Pl.l  Pi 

tl 

Pi , i^Pi^Pm  , i^Pi , i^Pi 

ft 


Pl  , 1 Pl^Pu,2  Pl^Pl.l  Pi 

ft 


^Pref  , » "Pdatum'^^Pref  , 2— Pdatum 

^Pref  .2^^'  ’ “Pdatum'^^Pref  . 1 -Pdatum 
^Pref  , • » "Pdatum'^^Pref  , 2— Pdatum 
^Pref  , 2^  * » "Pdatum'^^Pref  , 1— Pdatum 
^Pref  , 1""®  • » "Pdatum“^^Pref  , 2— Pdatum 
^Pref  , 2^®  • » "Pdatum'^^Pref  , 1— Pdatum 


139 

140 

141 

142 

143  * 

144  * 1 


yR=y2 

yL=y2 

yc==y2 


. 1 Pz^PjJ.Z  Pi 

n 


^L.  l^Pl=^U,2 

n 

Pi,i^Pi>Pv.z^Pz 

i 


^Pref 

1=0- 

» "Pdat 

^Pref 

2=0. 

» 'Pdat 

^Pref 

1-0- 

' "Pdat 

^Pref 

2=0- 

» "Pdat 

^Pref 

1=0. 

» “Pdat 

^Pref 

2-0. 

» "Pdat 

um'^^Pref  , 2— Pdatum 
um’^^Pref  , 1— Pdatum 
um^^Pref  , 2— Pdatum 
um“^^Pref  , 1— Pdatum 
um^^Pref  , 2— Pd  ATUM 
dm"^^Pref  , 1— Pdatum 
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TABLE  1:  Configurations /Environments  Used  to  Test  the  VENTCF  Subroutine 

(Cont ' d) 

= l.OOkg/m^  > P2  ""  0.50kg/m^  > pg  ==  0.25kg/m^ 
y,  = 0 . Om  < y,  = 1 . 5m  < y-  = 3 . Om  < y . 4 . 5m  < y,  =6  .Om 

= 1.0il25(10=)pa 


CASES  109-162:  SPACE  1 (TOP  SPACE)  IS  INSIDE  ROOM,  SPACE  2 (BOTTOM  SPACE) 

(CONT'D)  IS  OUTSIDE  SPACE  - APPENDIX:  FIGURE  Ic 


CASE  SPACE  SPACE 
NO.  1 2 


145 

146 

147 

148 

149 

150 


Upper  All 
Layer  One 

Layer 

yR=y3  yR=y4 
yL=y3  yL=y4 
yc=y5  yc“y4 


* 

* 


4. 


, 1 Pz^P-L  , 2 Pi 

n 

Pm.i^Pi^Pi^.z 

n 

P\}  . l^Pl^Ph  . Z^Pz 


^Pref.i  ® • » "Pdatum^^Pref , 2— Pdatum 
^Pref .z^^-  » "Pd  atum'^^Pref  i~Pdatum 
^Pref  , 1^®  • » "Pdatum'^^Pref  , 2— Pdatum 
^Pref  , 2""®  • » "Pdatdm^^Pref  . 1— Pdatum 
^Pref  , 1""®  * » "Pdatum'^^Pref  , 2— Pdatum 
^Pref  , 2^^  • » “Pdatum'^^Pref  , 1— Pdatum 


151 

152 

153 

154 

155  * 

156  * 


yR=y3 

yL=y3 

yc=y3 


i 


P\3  , i""P3^Pu . z""Pz^Pl  , z""Pi  ^Pref  , » “Pdatum'^^Pref  . 2— Pdatum 

^Pref  , 2""®  • » "Pdatum'^^Pref  , 1— Pdatum 
P\3 , i""Pz^P\j , z^Pi. . z""Pi  ^Pref  , • ’ "Pdatum'^^Pref  , 2— Pdatum 

^Pref  , 2^®  * ’ ’Pdatum^^Pref  , 1— Pdatum 
Pu,  , 2""P2^Pl  , 2""Pi  ^Pref,  i""®*  * "Pdatum'^^Pref  , 2-Pdatum 

^Pref  , 2^®  • > "Pdatum“^^Pref  , i-Pdatum 


157 

158 

159 

160 

161  * 
162  * i 


yR=y2 

yL==y2 

yc=y2 


Pu  . 1 Pz^P\3 . 2 Pi 

n 


Pm ,1  Pi  Pm .z 

tt 


^U . 1 ^l^Pu,2  Pz 

4. 


^Pref  , * » "Pdatum'^^Pref  , 2—PdATUM 

^Pref  , 2""®  • » "Pd atum^^Pref  , i-Pdatum 
^Pref  , • » "Pdatdm“^^Pref  , 2— Pdatum 

^Pref  , 2^^  • » "Pdatum'^^Pref  , 1— Pdatum 
^Pref  . » "Pdatum^^Pref  , 2-Pdatum 

^Pref  , 2^®  • » "Pdatum'^^Pref  , 1— Pdatum 


18 


Ap  (pascal) 


igure  1.  Plots  of  the  calculated  mass  flow  rate  to  upper  space  1,  i» 

and  to  lower  space  2,  functions  of  the  computed 

(coarse-scale)  cross-vent  pressure  difference,  Ap,  for  Cases  28, 
29,  and  30  (see  Table  1). 
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'^VENT,1 

(kg/s) 


'^VENT,2 

(kg/s) 


Figure  2.  Plots  of  the  calculated  mass  flow  rate  to  upper  space  1,  i, 

and  to  lower  space  2,  ^vent  2’  functions  of  the  computed  (fine- 
scale)  cross -vent  pressure  difference,  Ap,  for  Cases  28,  29,  and 
30  (see  Table  1) . 
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VENTCF  - 


APPENDIX 


Calculation  of  ti^e  Flow  Through  a Horizontal  Ceiling/Floor  Vent 
Connecting  Two  Spaces 
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VENTCF 


Calculation  of  the  Flow  Through  a Horizontal  Ceiling/Floor  Vent 
Connecting  Two  Spaces 


DESCRIPTION 


Consider  an  instant  of  time  during  the  simulation  of  a multi-room  compartment 
fire  environment.  This  algorithm  can  be  used  to  calculate  the  flow  of  mass, 
enthalpy,  oxygen,  and  other  products  of  combustion  through  a horizontal  vent 
located  in  a ceiling/floor  segment  common  to  any  two  inside  rooms  of  the 
facility  or  between  an  inside  room  and  the  outside  environment  local  to  the 
vent . 

Depicted  in  Figure  la  is  the  vent  and  the  two  spaces  when  they  are  both  inside 
rooms  of  a multi-room  facility.  Figures  lb  and  Ic  depict  the  situation  when 
the  two  spaces  involve  one  inside  room  of  the  facility  and  one  outside  space, 
either  above  or  below  the  room,  in  which  is  simulated  the  outside  environment 
local  to  the  vent. 

As  in  Figure  1,  designate  the  top  space  as  space  1 and  the  bottom  space  as 
space  2,  It  is  assumed  that  the  temperature,  density,  concentration  of  oxygen 
and  of  other  products  of  combustion  of  interest  in  the  upper  and  lower  layer 
of  each  inside  room  and  in  the  environment  local  to  the  vent  of  an  outside 
space  are  specified.  Also  specified  in  each  inside  room  are:  the  elevation 
above  a datum  elevation  of  the  floor,  and  the  upper- layer/lower- layer 
interface;  the  pressure  at  the  floor  of  each  room  above  a specified  datum 
pressure;  and  the  datum  absolute  pressure  there.  Specified  in  an  outside 
space  are:  a reference  elevation  above  the  datum  elevation,  and  the  pressure 
at  this  reference  elevation  above  the  specified  datum  pressure. 

The  calculation  of  flow  through  the  vent  is  based  on  the  model  of  reference 
[1]  which  takes  into  account  considerations  of  flow  due  to  a cross -vent 
pressure  difference  and,  when  appropriate,  of  flow  induced  through  the  vent 
when  existing  fluid  configurations  are  unstable,  i.e.,  when  a relatively  cool, 
dense  gas  in  the  upper  space  overlays  a less  dense  gas  in  the  lower  space. 
Whereas  a non-zero  cross -vent  pressure  difference  always  tends  to  drive  a 
unidirectional  flow  through  the  vent,  from  the  higher-  to  the  lower-pressure 
side,  flow  induced  by  an  unstable  fluid  density  configuration  tends  to  lead  to 
an  exchange  - type  of  flow,  with  fluid  in  the  lower  space  rising  into  the  upper 
space  and  flow  from  the  upper  space  dropping  into  the  lower  space.  As  a 
result,  the  flow  through  the  vent  can  be  unidirectional,  into  the  upper  or 
lower  space,  or  it  can  involve  a complex  two -component  flow,  whereby  flow  is 
drawn  simultaneously  from  the  bottom  to  the  top  space  and  from  the  top  to  the 
bottom  space. 

The  geometry  and  the  conditions  local  to  the  vent  which  determine  the 
characteristics  of  the  vent  flow  are  depicted  in  Figure  2.  These  include:  the 
densities,  and  p^,  and  the  hydrostatic  pressures,  p^  and  P2 , at  the 
elevation,  but  away  from  the  immediate  vicinity  of  the  vent  in  the  upper  and 
lower  spaces,  respectively,  and  the  area,  A^ , and  shape  of  the  vent. 

Regarding  the  shape,  at  the  present  time  results  for  horizontal  vent  flows  are 
only  available  for  circular  or  square  vents.  Other  properties  local  to  the 
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vent  and  indicated  in  Figure  2 are  and  T2  , the  absolute  temperatures,  Cq2  ^ 
and  Cq2^2’  concentrations  of  oxygen,  and  ^ and  Cj,  2 > ^ “ 2 to  Npj^Q^  , 

the  concentrations  of  a product  of  combustion  k’. 

In  the  present  algorithm,  the  calculation  of  the  vent  flow  is  based  on  the 
previously  developed  theory  presented  in  reference  [1] . When  the  scenario 
involves  an  unstable  density  configuration,  it  is  determined  there  that  an 
exchange  - type  of  flow  occurs  if  the  magnitude  of  the  cross -vent  pressure 
difference  is  less  than  that  of  a characteristic  value,  App^QQ^ , a specified 
function  of  cross -vent  density  difference  and  size  and  shape  of  the  vent. 

When  the  vent  flow  is  unidirectional  the  model  of  [1]  uses  the  incompressible 
orifice  - flow- type  of  flow  rate  calculation.  This  is  the  method  presented, 
e.g.,  in  reference  [2]  and  used  traditionally  in  most  zone-type  compartment 
fire  models  to  calculate  the  rate  of  flow  through  vents  of  vertical 
orientation.  In  order  to  be  able  to  treat  fire  scenarios  at  times  when  vent 
flows  are  driven  by  arbitrarily  high  cross-vent  pressure  differences,  i.e., 
when  compressibility  effects  begin  to  be  significant,  the  model  of  [1]  has 
been  modified  here  to  include  the  ideas  introduced  in  reference  [3]  and 
implemented  in  the  VENTHP  algorithm/subroutine  of  reference  [4] . Such  high 
cross -vent  pressure  differences  could  occur,  for  example,  in  fire  scenarios 
involving  flows  through  cracks  in  otherwise  hermetically- sealed  fire 
compartments . 

The  VENTCF  subroutine  has  been  subjected  to  extensive  parametric  testing  in  a 
"stand-alone"  mode.  This  is  described  in  reference  [5]. 


OUTPUT 


^02,1  t ^02 , 2 ] 

Concentration  of  oxygen  in  the  top  [bottom]  space  at 
the  elevation  of  the  vent,  [(kg  of  oxygen)/(kg  of 
layer) ] 

, 1 t , 2 ] 

Concentration  of  product  K in  the  top  [bottom]  space 
at  the  elevation  of  the  vent,  [(unit  of  product  K)/(kg 
of  layer) ] 


E N T , 0 2 , 


; I 


1 or  2 


Concentration  of  oxygen  in  the  component  of  the  vent 
flow  entering  space  I , provided  such  vent  flow 
component  is  non-zero,  i.e.,  provided  Invent,  !•  non- 
zero. [(kg  of  oxygen) /(kg  of  vent  flow)] 
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^VENT  K I’  1=1  or  2;  K=2  to  NPROD 

Concentration  of  product  of  combustion  K in  the 
component  of  the  vent  flow  entering  space  1 , provided 
such  vent  flow  component  is  nonzero,  i.e.,  provided 
^vENT.i-  non-zero,  [(unit  of  product)/(kg  of  vent 
flow)  ’] 


Mjj  I [Ml  I ] ! ^ ^ 1 2 

If  space  1 is  an  inside  room:  Rate  at  which  mass  is 
added  to  the  upper  [lower]  layer  of  room  1 due  to  the 
vent  flow  component  which  enters  the  other  space. 

Note  that  this  will  alv  s be  negative,  since  the 
layer  which  supplies  tht  material  flowing  to  the  other 
space  will  always  have  mass  extracted  from  it.  [W] 

If  space  1 is  an  ’’outside  room”  : M^  j is  the  rate  at 
which  mass  is  added  to  space  1 due  to  the  vent  flow 
component  which  enters  the  other  space  (i.e.,  the 
inside  room) . Note  that  this  will  always  be  negative, 
since  an  outside  space  which  supplies  material  flowing 
through  the  vent  to  an  adjacent  space  will  always  have 
mass  extracted  from  it.  M^  j is  identical  to  My  j . 

[W] 


M • 1 = 1 to  2 

^•^VENT  I’  ^ 

Mass  flow  rate  of  the  vent  flow  component  entering 
space  1.  [kg/s] 


^0  2 , L , I t ^0  2 , U , I 


1 or  2 


If  space  1 is  an  inside  room:  Rate  at  which  oxygen  is 
added  to  the  upper  [lower]  layer  of  room  1 due  to  the 
vent  flow  component  which  enters  the  other  space. 

Note  that  this  will  always  be  negative,  since  the 
layer  which  supplies  the  material  flowing  through  the 
vent  to  the  other  space  will  always  have  its  convected 
oxygen  extracted  from  it.  [(kg  of  02)/s] 

If  space  1 is  an  ’’outside  room":  Pq2  y j is  the  rate 
at  which  oxygen  is  added  to  space  1 due  to  the  vent 
flow  component  which  enters  the  other  space  (i.e.,  the 
inside  room) . Note  that  this  will  always  be  negative, 
since  an  outside  space  which  supplies  material  flowing 
through  the  vent  to  an  adjacent  space  will  always  have 
its  convected  oxygen  extracted  from  it.  Pq2  l i 
identical  to  Pq2  u I'  [ O2VS] 
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1 to  2 


^0  2 , V E N T , I > 

Mass  flow  rate  of  oxygen  of  the  vent  flow  component 
entering  space  I.  [(kg  of  oxygen)/s] 

^K,L,i  ^ 2 to  I = 1 or  2 

If  space  I is  an  inside  room:  Rate  at  which  product 
of  combustion  K is  added  to  the  upper  [lower]  layer  of 
room  I due  to  the  vent  flow  component  which  enters  the 
other  space.  Note  that  this  will  always  be  negative 
since  the  layer  which  supplies  the  material  flowing 
through  the  vent  to  the  other  space  will  always  have 
its  convected  product  of  combustion  K extracted  from 
it.  [(unit  of  product  K)/s] 

If  space  I is  an  "outside  room":  u i rate 

at  which  product  of  combustion  K is  added  to  space  I 
due  to  the  vent  flow  component  which  enters  the  other 
space  (i.e.,  the  inside  room).  Note  that  this  will 
always  be  negative,  since  an  outside  space  which 
supplies  material  flowing  through  the  vent  to  an 
adjacent  space  will  always  have  the  convected  product 
of  combustion  K extracted  from  it.  Pj,  ^ i is 
identical  to  ^ j.  [(unit  of  product  K)/s] 

^K,  VENT,  I’  i i 2 to  NppQjj 

Flow  rate  of  product  of  combustion  K in  the  vent  flow 
component  entering  space  I.  [(unit  of  product  K)/s] 

Pi  [Pz] 

Absolute  hydrostatic  pressure  in  the  top  [bottom] 
space  at  the  elevation  of  the  vent.  [Pa  = kg/(m-s^)] 

Qu,i  I = 1 or  2 

If  space  I is  an  inside  room:  Rate  at  which  enthalpy 
added  to  the  upper  [lower]  layer  of  room  I due  to  the 
vent  flow  component  which  enters  the  other  space. 

Note  that  this  will  always  be  negative,  since  the 
layer  which  supplies  the  material  which  flows  to  the 
other  space  will  always  have  its  convected  enthalpy 
extracted  from  it.  The  enthalpy  is  based  on  the 
absolute  temperature  of  the  flow,  Tyj.j,.j.  j.  [W] 
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If  space  I is  an  "outside  room" : Qy  j is  the  rate  at 
which  er  r.alpy  is  added  to  space  I due  to  the  vent 
flow  component  which  enters  the  other  space  (i.e,,  the 
inside  room')  . Note  that  this  will  always  be  negative, 
since  an  o side  space  which  supplies  material  which 
flows  througn  the  vent  to  an  adjacent  space  will 
always  have  its  convected  enthalpy  extracted  from  it. 
The  enthalpy  is  based  on  the  absolute  temperature  of 
the  flow,  Ty  ^ . Qt  j is  iaentical  to  Qy  j . [W] 


Qvznt, I > I = 1 to  2 

Total  enthalpy  flow  rate  in  the  vent  flow  component 
entering  space  I.  This  is  based  on  the  absolute 
temperature  of  the  flow,  j . [W] 

Ti  [Tj] 

Absolute  temperature  in  the  top  [bottom]  space  at  the 
elevation  of  the  vent.  [K] 


^VENT.I’  ^ lto2 

Absolute  temperature  of  the  vent  flow  component 
entering  space  I , provided  such  vent  flow  component  is 
non-zero.  [K] 

Ap 

P2  "Pi.  i.e.,  pressure  in  the  bottom  space  at  the 
elevation  of  the  vent  - pressure  in  the  top  space  at 
the  elevation  of  the  vent.  [Pa  - kg/(m-s^)] 

Pi  [Pi] 

Density  in  the  top  [bottom]  space  at  the  elevation  of 
the  vent.  [kg/m^’^S] 


PvENT  , I ’ ^ ^ ^ 


Density  of  the  vent  flow  component  entering  space  I, 
provided  such  vent  flow  component  is  nonzero,  [kg/m^ ] 
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INPUT 


Av 


Area  of  the  vent  [m^ ] . 


, K , I ^ 2 to  Npj^Qp  ; I 1 or  2 

If  space  I is  an  inside  room:  Concentration  of 
product  of  combustion  K in  lower  [upper]  layer  of  room 
I if  the  volume  of  the  lower  [upper]  layer  is  non-zero 
(if  the  lower  [upper]  layer  volume  is  zero  then  the 
K I K I ^ value  is  not  used  in  the  calculation) . 

[(unit  of  product  K)/(kg  of  layer)] 

If  space  I is  an  "outside  room":  j is  the 

uniform  concentration  throughout  the  space  of  product 
of  combustion  K;  c^  j is  specified  as  being 
identical  to  c^  j.  j . [(unit  of  product  K)/(kg  of  local 
atmosphere) ] 


, 02  , I 


'U  , 02  , I 


1 or  2 


If  space  I is  an  inside  room:  Concentration  of  oxygen 
in  lower  [upper]  layer  of  room  I if  the  volume  of  the 
lower  [upper]  layer  is  non- zero  (if  the  lower  [upper] 
layer  volume  is  zero  then  the  q2  i 02  i ^ value 

is  not  used  in  the  calculation),  [(leg  of  oxygen) /(kg 
of  layer) ] 

If  space  I is  an  "outside  room":  c^  02  i 
uniform  concentration  throughout  the  space  of  oxygen; 
Cu  02  I specified  as  being  identical  to  c^  02  i • 
[(kg  of  oxygen)/(kg  of  local  atmosphere)] 


Specific  heat  at  constant  pressure  of  the  vent  flow. 
[W-s/(kg'K)]  (suggest  10^  W's/(kg-K)  for  air  as 
default) 


Characteristic  dimension  of  the  vent:  the  diameter  of 
a circular  vent  and  the  length  of  a side  of  a square 
vent,  [m] 
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Maximum  allowed  number  of  products  of  combustion. 


Np 


MAX 


N 

PROD 


Number  of  products  of  combustion,  including  oxygen, 
being  tracked  in  the  simulation. 


N 


s 


RAPE 


Number  associated  with  the  shape  of  the  horizontal 
vent  (the  present  algorithm  is  based  on  results  for 
circular  and  square  vents) 


N 


SHAPE 


1 for  a circular  vent 

2 for  a square  vent 


Pdatum 


Datum  absolute  pressure.  [Pa  = kg/(m-s^)] 


, I t ^ I ] ; 1 1 or  2 

If  space  I is  an  inside  room:  Absolute  temperature  of 
the  lower  [upper]  layer  in  room  I if  the  volume  of  the 
lower  [upper]  layer  is  non- zero  (if  the  lower  [upper] 
layer  volume  is  zero  then  the  j i 1 is  not 

used  in  the  calculation) . [K] 

If  space  I is  an  ’’outside  room”:  j is  the  uniform 

absolute  temperature  there,  taken  to  be  the 
temperature  at  the  reference  elevation,  yj^^p  i ! *^0  i 
is  specified  as  being  identical  to  j . [K] 

yCEIL.I  ^yREF,I^ 


If  space  I is  an  inside  room:  Elevation  of  the 
ceiling  [floor]  of  room  I above  the  datum  elevation, 
[m] 

If  space  I is  an  "outside  room":  Yceil,!  ^ref,! 

are  both  identical  and  equal  to  the  reference 
elevation  of  space  I above  the  datum  elevation,  i.e., 
the  specification  must  satisfy  y^piL  i = Yref,!- 
latter  identity,  which  will  never  be  satisfied  for  an 
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A YER , I ’ 


yVEN  T 


•^Pref  , I ’ ^ 


inside  room,  is  a characteristic  of  the  input  data 
used  to  distinguish  an  inside  room  from  an  "outside 
room . " [m] 


= 1 or  2 

If  space  I is  an  inside  room:  Elevation  of  the 
upper/lower  layer  interface  in  room  I above  the  datum 
elevation,  [m] 

If  space  I is  an  "outside  room":  Yl^yer  i specified 
as  being  identical  to  j . [m] 


Elevation  of  the  vent  above  the  datum  elevation.  Note 
that  Yvent  be  identical  to  either  Yceil  i 

Yref  I for  each  of  the  one  or  two  inside  rooms 
involved  in  the  calculation,  [m] 


= 1 or  2 

Pressure  at  the  reference  elevation,  Yref  i » space 
I above  the  datum  absolute  pressure,  Pd^tum • space 

I is  an  inside  room,  then  5pj^pp  j and  Yref  i ^^^t 
correspond  to  the  pressure  and  elevation, 
respectively,  within  the  room  and  at  the  floor.  [Pa  = 
kg/(m- s^ ) ] 


Error  tolerance  for  ^Pref  i * Perror  i defined 

as  the  uncertainty  in  ^Pj^^p  j,  I =■  1 or  2 , then 
Perror, I satisfies 

IPerror.iI  ^ ^pPo  **”  I ^Pref  * I I 

where  Pq  = 1.0  Pa.  The  first  term  is  based  on  an 
absolute  error  tolerance  and  dominates  the  above  error 
bound  when  |5Pj^ef  i I less  than  1.0  Pa.  The  second 
term  is  a relative  error  tolerance  and  dominates  when 
I^Pref  I I greater  than  1.0  Pa.  should  be  chosen 

to  be  consistent  with  the  tolerance  specified  for  the 
computation  of  5pj^pp  j terms  in  the  overall 
compartment  fire  model  computer  code  which  uses  this 
algorithm. 
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Pl , I t , I ] ’ ^ ^ ^ 

If  space  I is  an  inside  room:  Density  of  the  lower 
[upper]  layer  in  room  i if  the  volume  of  the  lower 
[upper]  layer  is  non-zero  (if  the  lower  [upper]  layer 
volume  is  zero  then  the  Pl  i i ^ value  is  not  used 

in  the  calculation),  [kg/m^ ] 

If  space  I is  an  "outside  room" : i is  the  uniform 

density  there;  j is  specified  as  being  identical  to 
Pl  I-  [kg/m^] 

CALCULATION 

1.  Calculate  pj  for  1=1  and  2 , and  then  Ap  (the  pj  calculation  follows 
the  DELP  algorithm/subroutine  of  reference  [6]): 


Pi  ^Pref.i  ^Pi  Pdatum 

Ap  = P2  ■ Pi  ^ ^^PrEF,2  ’ ^PrEF,1^  ^^Pz  ” *^£1^ 

where 

Nref , I 

5pi 

else : 

5Pi  = ■ Pl , I sCYlayer , I ■ Yref.i) 

(3) 


- Yvent  - Ylayer,!  Pken: 

^ ■ Pl,iS(Yvent  ■ Yref.i^ 


and  where  g,  the  acceleration  of  gravity,  is  9.8  m/s^ . 

2 . Calculate  pj  , , ^02  i * and  the  Cjr  j , K = 2 to  ^prod  * for  I = _ iid  2 , 

and  then  Ap : 
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if  { [ (Yvent 


Yref.i)  ^Ylayer,!  Yref.i^^ 


[ ^Yvent  ycEiL.i^  ‘^Ylayer,!  YcEiL.i^]^  them 

Pi  , I ’ f U , I » ^^0  2,1  ^U,0  2,I’  , I ^ ^U  , K , I 

else : 

Pi  ^ ^L  , I ’ ""  f L , I ’ ^02,1  ^ ^L,02,I>  , I ""  ^L  , K . I 

(4) 

4p  = Pi  - p2  (5) 

3.  Define  Vg^  j as  the  volume  rate  of  flow  through  the  vent  into  space  I 
that  is  predicted  with  a "standard,"  unidirectional -flow- type 
calculation,  where  arbitrarily  high  cross-vent  pressures  are  allowed. 
Here,  the  calculation  follows  the  model  of  reference  [3]  as  implemented 
in  the  VENTHP  algorithm/subroutine  of  reference  [4] . Calculate  Vg^  j , I 
= 1 and  2 : 


If  Ap  = 0 then; 


V 


S T , 1 


= V 


S T , 2 


= 0 


(6) 


If  Ap  > 0 or  Ap  < 0 then: 

Define  and  compute  p,  e,  and  x; 

If  Ap  > 0,  i.e.,  "standard"  flow  from  (lower)  space  2 to 
(upper)  space  1,  then: 

P = P2',  € = Ap/p2  (7) 

If  Ap  < 0,  i.e.,  "standard"  flow  from  (upper)  space  1 to 


(lower)  space  2,  then: 

p = p^\  € = -Ap/pi 

(8) 

X = 1 - € 

(9) 

Compute  C(x),  the  vent  flow  coeffient,  and  w(x) : 
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w(x) 


(11) 


r 1 - [3/(47)]€  if  0 < e < 10"^ 

1 f(x)/[2£]i/2  if  1 > £ > 10-5 


where 


f(x) 


{[27/(7  - - x^^  ■ l)/7j)l/2 

if  € < 1 - [2/(7  + 

{7[2/(7  + 1)]^^  + l)/(7  - 1) )l/2 

if  € > 1 - [2/(7  + 1)]^/ 


1 ) 

(12) 

1 ) 


and  where  7,  the  ratio  of  specific  heats  of  the  vent  flow  gas,  is 
taken  to  be  that  of  air,  1.40. 

Define  and  compute  Ap|i^^  : 

^Pc<JT  [€p pascal,  l<5pj^j.p  J , I 5pj^j.p  2 1 ) ]^^^  (13) 

Define  and  compute  F^oise’  ^ numerical  damping  factor: 

^NOiSE  = 1-0  - axp(- |Ap|^/ VApJ^;?)  (14) 


Compute  V,  the  "standard"  volume  flow  rate  from  the  high  to  low 
pressure  space,  and  j,  1=1  and  2: 

^ = Fj,0PSpC(x)w(x)(2/p)'/2^vl^p|i/2  (15) 

VsT,i  = V,  VsT,2  = 0 if  Ap  > 0 (16) 

^ST,2  = "^ST,!  = ^ if  Ap  < 0 (17) 

Note  in  Eqs . (7)  and  (8)  that  e -»•  0 , i.e.,  x -*•  1 , corresponds  to 

the  incompressible  limit  where  C -*■  0.68  in  Eq.  (10)  according  to 
the  recommendation  of  reference  [2]. 
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The  term  Fj^Qjgg  of  Eq . (14)  is  designed  to  damp  out  the  numerical 
noise  (error)  in  the  calculated  value  for  Ap  that  would  otherwise 
be  dominant  in  Eq . (15)  when  Ap  is  small  relative  to  the  maximum 
of  1.0  pascal  and  the  calculated  reference  pressures,  SpR^p  i, 
5prep  2-  term  Ap^^^^  of  Eq.  (13)  is  an  estimate  of  how  small 

the  maximum  of  | Ap | must  be  to  retain  a few  digits  of  accuracy  in 
the  calculation  of  Ap . When  the  calculated  value  of  |Ap|  is 
smaller  than  P^ut > value  and,  therefore,  the  value  of  V in 

Eq . (15)  will  likely  contain  noise  which  should  be  damped. 

is  constructed  to  tend  towards  1 when  | Ap | is  large  relative  to 
ApcuT  tends  towards  0 when  1 Ap  1 is  small  relative  to  Ap^^^  . 


Define  and  calculate  as  the  volume  rate  of  flow  exchanged  through 

the  vent  between  the  two  spaces  because  of  an  unstable  cross -vent 
density  configuration,  i.e.,  when  Ap  > 0: 


If  Ap  < 0:  = 0 


(18) 


If  Ap  > 0: 


Calculate  AppLooD • 

^PfLOOD  APE  ) 

where 


0.754  for  circular  vents 


'SHAPE 


•« 


0.942  for  square  vents 

Calculate  , the  maximum  value  of  : 

^EX.MAX  = 0.10[2gApAv5'V(Pi  + 
Calculate  : 


(19) 


(20) 


(21) 


^EX.MAX^^  ' I 1 / I ^PpLOOD I ^ 
0 Otherwise 


if  1 Ap  I / I AppEooD  I ^ 

(22) 
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5.  Define  and  calculate  j,  I = 1 and  2,  as  the  volume  rate  of  the 

vent  flow  component  entering  space  I . 


E N T , I "^S  T , I "'^E  X 


(23) 


6.  Calculate  the  vent  flow  properties  Pvent,i>  ^vent,02,i- 


'VENT  , K , I 


, K=2  toNpj^Qp,  1 = 1 and  2: 


^VENT  ,1  ^2  > ^VENT  , 2 


T = T T = T 

-^VENT  ,1  -^2  ’ -^VENT  ,2 


(24) 

(25) 


ENT,  02,1  ^0  2,2’  ENT,  02,  2 ^0  2,1 


(26) 


^VENT,K,l  ^K,2’  ^VENT,K,2  ^K,l’  ^ ^^°^PR0D 


(27) 


Calculate  the  vent  flow  rates  M. 


•V  E N T , I 


,1=1  and  2 : 


ENT,  I ^VENT.l'^VENT,! 


(28) 


Calculate  the  vent  flow  rates  Qvent,!’  ^02,vent,I’  ^k,vent,i,  K = 2 to 

^PROD  ’ ^ ^ ^ 


Qvent,  I ^vent  , i ‘^p'^vent  , i 


^02, vent, I ^VENT,I^VENT,02,I 


(29) 

(30) 


^K,VENT,^  “ ^VENT,I^VENT,K,I’  ^ ^ ^PROD 


(31) 


Calculate  rates  at  which  flows  are  added  to  layers  of  each  space  as  a 
result  of  the  vent  flow  extracted  from  it.  First  consider  space  1 (I  = 
1,  J = 2)  and  then  space  2 (1=2,  J = 1) . For  either  case: 
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HCYveNT  “ yREF.I^  ^yLAYER.I  “ yREF.I^^ 

t^yVENT  ""  yCEIL.l)  (yLAYER,!  ^ yCEIL.I^]^ 

(i.e.,  the  vent  flow  to  room  J is  extracted  from  the  upper  layer  of  room 
I and,  if  space  I is  an  inside  room,  the  lower  layer  of  room  I is 
unchanged)  then: 


M, 


u , I 


02  , U , I 


K , U , I 


If  yR] 


- M 

^^VENT  , J > 

Ml.i  - 

0 

(32) 

Qv  E N T , J ’ 

Ql.i  = 

0 

(33) 

fo2,VENT,J’ 

^0  2 , L , I ® 

(34) 

, V E N T , 

J ’ 

O 

n 

M 

(35) 

, I yCE  I L , I 

the  results 

(i.e. 
of  Eq . 

, space  I is  an  outside 
(32)-(35)  as  follows: 

space) 

t.I  = ^U.I 

(36) 

. I = Qu  . I 

(37) 

' = p 

02 , L , I ^02 

, u . I 

(38) 

K , L , I , U 

, I 

(39) 

If  the  condition  above  Eq.  (32)  is  not  satisfied 

(i.e.,  the  vent  flow  to  room  J is  extracted  from  the  lower  layer  of  room 
I and,  if  space  I is  an  inside  room,  the  upper  layer  of  room  I is 
unchanged)  then: 


^L.I  = 

- M 

^^VENT  , J » 

Mu, I =0 

(40) 

!l,i  = 

Qv  E N T . J ’ 

Qu , I - 0 

(41) 

*0  2 , L , I 

= - P P =0 

^02  , VENT  , J » ^02  , U , I ^ 

(42) 

K , L , I 

. V E N T ^ 

P =0 

, J » , U , I ^ 

(43) 

If  yREF  I ycEiL  I space  I is  an  outside  space) 

modify  the  results  of  Eqs . (40) -(43)  as  follows: 
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SUBROUTINE  VARIABLES 


All  nomenclature  in  the  subroutine  is  identical  to  the  nomenclature  used  above 
except  for: 

Av  - AVENT  [m^] 
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c 

COEF 

s 

CP  [W-s/(kg-K)] 

r 

'^SHAPE 

CSHAPE  [dimensionless] 

, K , I > , K . I 

CONL(K,I),  CONU(K,I),  I = 1 or  2 [(unit  of 
product  K)/(kg  of  layer)] 

, 02 , I ’ , 02 , I 

C0NL(1,I),  C0NU(1,I),  I - 1 or  2 [(kg  of 
oxygen) /(kg  of  layer)] 

^VENT  , K , I 

CVENT(K,I),  1=  1 or  2 [(unit  of 
product)/(kg  of  vent  flow) 

^VENT  , 02  , I 

CVENT(1,I),  I = 1 or  2 [(kg  of  oxygen) /-(kg 
of  vent  flow) ] 

, I 

C(K,I),  1=  1 or  2 [(unit  of  product)/(kg 
of  vent  flow) 

^02,1 

C(1,I),  I - 1 or  2 [(kg  of  oxygen)/(kg  of 
vent  flow) ] 

0 I S E 

FNOISE  [dimensionless] 

f 

FF  [dimensionless] 

g 

9 . 8 m/s^ 

M M 

, I ’ , I 

XML(I),  XMU(I),  I - 1 or  2 [kg/s] 

M 

“vent  , I 

XMVENT(I)  [kg/s] 

N 

^^PMAX 

NPMAX  [dimensionless] 

N 

PROD 

NPROD  [dimensionless] 

N 

“shape 

NSHAPE  [dimensionless] 

, L , I ’ , U , I 

PL(K,I),  PU(K,I),  I - 1 or  2 [(unit  of 
product  K)/s] 

, V E N T , I 

PVENT(K,I),  I - 1 or  2 [(unit  of  product 
K)/s] 

P P 

^02  . L , I ’ ^02  , U , I 

PL(1,I),  PU(1,I),  I - 1 or  2 [(kg  of 
oxygen) /s] 

^0  2 , V E N T , I 

PVENT(1,I),  I = 1 or  2 [ (kg  of  oxygen)/s] 

Pd  ATUM 

PDATUM  [Pa  = kg/(m*s2)] 
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Pi 

- 

P(I) , I = 1 or  2 [Pa  = kg/(m-s^)] 

Ql , I ' Qu , I 

- 

QL(I) , QU(I) , I = 1 or  2 [W] 

Qv  E N T , I 

- 

QVENT(I),  I = 1 or  2 [W] 

Ti 

- 

T(I) , I = 1 or  2 [K] 

'^L  . I ’ , I 

- 

TL(I) , TU(I) , I = 1 or  2 [K] 

E N T , I 

- 

TVENT(I) , I = 1 or  2 [K] 

V 

- 

V [m^/s] 

^EX 

- 

VEX  [mVs] 

^E  X , M A X 

- 

VEXMAX  [mVs] 

T , I 

- 

VST(I)  , I - 1 or  2 [mVs] 

E N T , I 

- 

WENT(I)  , I - 1 or  2 [mVs] 

W 

- 

W [dimensionless] 

X 

- 

X [dimensionless] 

yREF  , I 

- 

YREF(I) , I = 1 or  2 [m] 

Flayer , i 

- 

YLAY(I) , I = 1 or  2 [m] 

yvENT 

- 

YVENT  [m] 

yCE I L , I 

- 

YCEIL(I),  I = 1 or  2 [m] dimensionless ] 

7 

- 

1.40 

Ap 

- 

DEEP  [Pa  = kg/(m*s2)] 

^Pc<JT 

- 

DPC1D2  [Pa  - kg/(m-s2)] 

^Pflood 

- 

DELPFD  [Pa  = kg/Cm-s^)] 

Ap 

- 

DELDEN  [kg/m^] 

5pi 

- 

DP(I) , I = 1 or  2 [Pa  = kg/(m-s^)] 

"^Pref  , I 

- 

DPREF(I),  I = 1 or  2 [Pa  = kg/Cm-s^)] 

€ 

- 

EPS  [dimensionless] 

- 

EPSP  [dimensionless] 
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p - RHO  [kg/m^] 

Pj  - DEN(I),  I = 1 or  2 [kg/m^ ] 

Pu.i  ■ DENL(I),  DENU(I),  I = 1 or  2 [kgV  ] 

Pvent.i  - DENVNT(I),  I = 1 or  2 [kg/m^ ] 

PREPARED  BY 

Leonard  Y.  Cooper 
October  1990 
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Figure  1.  The  possible  configurations  of  the  two  spaces  joined  by  a 

horizontal  ceiling/floor  vent  with  space  1 above  space  2:  a)  two 
inside  rooms;  b)  an  outside  space  above  an  inside  room;  c)  an 
inside  room  over  an  outside  space. 
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1 


Figure  2. 


The  geometry  and  conditions  local  to  a horizontal  ceiling/floor 
vent  which  determine  the  characteristics  of  the  vent  flow. 
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SUBROUTINE  VENTCF 


SUBROUTINE  VENTCF ( 

I AVENT , CONL , CONU , CP , D , NPMAX , NPROD , NSHAPE , PDATUM , TL , TU , 

I YCEIL,YREF,YLAY,YVENT,DPREF,EPSP,DENL,DENU, 

0 C,CVENT,XML,XMU,XMVENT,PL,PU,PVENT,P,QL,QU,QVENT,T,TVENT, 

0 DEEP, DEN, DENVNT) 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

C*BEG 

C***  VENTCF  - CALCULATION  OF  THE  FLOW  OF  MASS,  ENTHALPY,  OXYGEN  AND 
C OTHER  PRODUCTS  OF  COMBUSTION  THROUGH  A HORIZONTAL 

VENT  JOINING  AN  UPPER  SPACE  1 TO  A LOWER  SPACE  2.  THE 
SUBROUTINE  USES  INPUT  DATA  DESCRIBING  THE  TWO -LAYER 
ENVIRONMENT  OF  INSIDE  ROOMS  AND  Tl.E  UNIFORM 
ENVIRONMENT  IN  OUTSIDE  SPACES. 

SUBROUTINE  ARGUMENTS 


INPUT 


AVENT 

CONL(K,I) 


CONL(K,I) 


CP 
D 

NPMAX 
NPROD 
NSHAPE 

PDATUM 
TL(I) 

TU(I) 

YCEIL(I) 

YREF(I) 

YLAY(I) 

C YVENT 
C DPREF(I) 
C 

C EPSP 
C DENL(I) 

C 

C DENU(I) 


AREA  OF  THE  VENT  [M**2] 

CONCENTRATION  OF  PRODUCT  K IN  LOWER  LAYER  OF  AN  INSIDE 
ROOM  I OR  UNIFORM  CONCENTRATION  OF  PRODUCT  K IN  AN 
OUTSIDE  SPACE  I [(UNIT  OF  PRODUCT)/(KG  LAYER)] 
CONCENTRATION  OF  PRODUCT  K IN  UPPER  LAYER  OF  AN  INSIDE 
ROOM  I OR  UNIFORM  CONCENTRATION  OF  PRODUCT  K IN  AN 
OUTSIDE  SPACE  I [(UNIT  OF  PRODUCT)/(KG  LAYER)] 

SPECIFIC  HEAT  [W*S/(KG*K) ] 

DIAMETER  OF  VENT  FOR  CIRCLE,  LENGTH  OF  SIDE  FOR  SQUARE 
[M] 

MAXIMUM  ALLOWED  NUMBER  OF  PRODUCTS 

NUMBER  OF  PRODUCTS  IN  CURRENT  SCENARIO 

NUMBER  CHARACTERIZING  VENT  SHAPE:  1 = CIRCLE,  2 = 

SQUARE 

DATUM  ABSOLUTE  PRESSURE  [PA  = KG/(M*S**2)] 

TEMPERATURE  OF  LOWER  LAYER  OF  AN  INSIDE  ROOM  I OR 
TEMPERATURE  OF  AN  OUTSIDE  SPACE  I [K] 

TEMPERATURE  OF  UPPER  LAYER  OF  AN  INSIDE  ROOM  I OR 
TEMPERATURE  OF  AN  OUTSIDE  SPACE  I [K] 

HEIGHT  OF  CEILING  ABOVE  DATUM  ELEVATION  FOR  AN  INSIDE 
ROOM  I OR  YREF(I)  FOR  AN  OUTSIDE  SPACE  I [M] 

HEIGHT  OF  REFERENCE  ELEVATION  FOR  SPACE  I ABOVE  DATUM 
ELEVATION  [M] 

HEIGHT  OF  LAYER  ABOVE  DATUM  ELEVATION  FOR  AN  INSIDE 
ROOM  I OR  YREF(I)  FOR  AN  OUTSIDE  SPACE  I [M] 

HEIGHT  OF  VENT  ABOVE  DATUM  ELEVATION  [M] 

PRESSURE  IN  SPACE  I AT  ITS  REFERENCE  ELEVATION  ABOVE 
DATUM  ABSOLUTE  PRESSURE  [PA  = KG/(M*S**2)] 

ERROR  TOLERANCE  FOR  DPREF  [DIMENSIONLESS] 

DENSITY  OF  LOWER  LAYER  OF  AN  INSIDE  ROOM  I OR  DENSITY 
OF  AN  OUTSIDE  SPACE  I [KG/M**3] 

DENSITY  OF  UPPER  LAYER  OF  AN  INSIDE  ROOM  I OR  DENSITY 
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OF  AN  OUTSIDE  SPACE  I [KG/M**3] 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


OUTPUT 

C(K,I) 

CVENT(I) 

XML(I) 

XMU(I) 

XMVENT(I) 

PL(K,I) 

PU(K,I) 

PVENT(K,I) 

P(I) 

QL(I) 

QU(I) 

QVENT(I) 

T(I) 

TVENT(I) 

DELP 

DEN(I) 

DENVNT(I) 


- CONCENTRATION  OF  EACH  PRODUCT  IMMEDIATELY  ABOVE  (IN 
SPACE  1=1)  AND  BELOW  (IN  SPACE  1=2)  THE  VENT 
[(UNIT  OF  PRODUCT) /(KG  LAYER)] 

- CONCENTRATION  OF  EACH  PRODUCT  IN  THE  VENT  FLOW 
COMPONENT  ENTERING  SPACE  I [(UNIT  OF  PRODUCT) /(KG  OF 
VENT  FLOW)] 

- RATE  AT  WHICH  MASS  IS  ADDED  TO  THE  LOWER  LAYER  OF  AN 
INSIDE  ROOM  1=1  (1=2)  OR  TO  AN  OUTSIDE  SPACE  I = 

1 (I  = 2)  DUE  TO  THE  VENT  FLOW  COMPONENT  ENTERING 
SPACE  = 2 (I  = 1)  [KG/S] 

- RATE  AT  WHICH  MASS  IS  ADDED  TO  THE  UPPER  LAYER  OF  AN 
INSIDE  ROOM  I = 1 (I  = 2)  OR  TO  AN  OUTSIDE  SPACE  1=1 
(1=2)  DUE  TO  THE  VENT  FLOW  COMPONENT  ENTERING  SPACE 
I = 2 (I  = 1)  [KG/S] 

- MASS  FLOW  RATE  IN  THE  VENT  FLOW  COMPONENT  ENTERING 
SPACE  I [KG/S] 

- RATE  AT  WHICH  PRODUCT  K IS  ADDED  TO  THE  LOWER  LAYER  OF 
AN  INSIDE  ROOM  I = 1 (I  = 2)  OR  TO  AN  OUTSIDE  SPACE  I 
= 1 (I  = 2)  DUE  TO  THE  VENT  FLOW  CONPONENT  ENTERING 
SPACE  I = 2 (I  = 1)  [(UNIT  OF  PRODUCT)/S] 

- RATE  AT  WHICH  PRODUCT  K IS  ADDED  TO  THE  UPPER  LAYER  OF 
AN  INSIDE  ROOM  I = 1 (I  = 2)  OR  TO  AN  OUTSIDE  SPACE  I 
= 1 (I  = 2)  DUE  TO  THE  VENT  FLOW  COMPONENT  ENTERING 
SPACE  1=2  (1=1)  [(UNIT  OF  PRODUCT)/S] 

- FLOW  RATE  OF  PRODUCT  K IN  THE  VENT  FLOW  COMPONENT 
ENTERING  SPACE  I [(UNIT  OF  PRODUCT) /S] 

- ABSOLUTE  PRESSURE  IMMEDIATELY  ABOVE  (IN  SPACE  1=1) 
AND  BELOW  (IN  SPACE  1=2)  THE  VENT  [PA  = KG/(M*S**2)] 

- RATE  AT  WHICH  ENTHALPY  IS  ADDED  TO  THE  LOWER 

LAYER  OF  AN  INSIDE  ROOM  I = 1 (I  = 2)  OR  TO  AN  OUTSIDE 
SPACE  I = 1 (I  = 2)  DUE  TO  THE  VENT  FLOW  COMPONENT 
ENTERING  SPACE  I = 2 (I  = 1)  [W] 

- RATE  AT  WHICH  ENTHALPY  IS  ADDED  TO  THE  UPPER  LAYER  OF 
AN  INSIDE  ROOM  I = 1 (I  = 2)  OR  TO  AN  OUTSIDE  SPACE  I 
= 1 (I  = 2)  DUE  TO  THE  VENT  FLOW  COMPONENT  ENTERING 
SPACE  I = 2 (I  = 1)  [W] 

- FLOW  RATE  OF  ENTHALPY  IN  THE  VENT  FLOW  COMPONENT 
ENTERING  SPACE  I [W] 

- ABSOLUTE  TEMPERATURE  IMMEDIATELY  ABOVE  (IN  SPACE  I = 

1)  AND  BELOW  (IN  SPACE  1=2)  THE  VENT  [K] 

- ABSOLUTE  TEMPERATURE  OF  THE  VENT  FLOW  COMPONENT 
ENTERING  SPACE  I [K] 

- CROSS-VENT  PRESSURE  DIFFERENCE,  P(2)  - P(l) , [PA  = 
KG/(M*S**2)] 

- DENSITY  IMMEDIATELY  ABOVE  (IN  SPACE  1=1)  AND  BELOW 
(IN  SPACE  1=2)  THE  VENT  [KG/M**3] 

- DENSITY  OF  THE  VENT  FLOW  COMPONENT  ENTERING  SPACE  I 
[KG/M**3] 
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c 

C*END 

DIMENSION  CONL(NPMAX,*) , CONU(KPMAX,*) , TL(*) , TU(*) 

DIMENSION  YCEIL(*) , YREF(*) , YIAY(*) , DPREF(*) 

DIMENSION  DENL(*) , DENU(*) 

DIMENSION  C(NPMAX,*),  CVENT(NPMAX, *) , XML(*) , XMU(*) 

DIMENSION  XMVENTC*) 

DIMENSION  PL(NPMAX,*),  PU(NPMAX,*),  PVENT(NPMAX, *) , P(*) 

DIMENSION  QL(*),QU(*) 

DIMENSION  QVENT(*),  T(*) , TVENT(*) , DEN(*) , DENVNT(*) 

DIMENSION  DP(2),  VST(2)  , WENT(2) 

PARAMETER  (GAM=1.40D0) 

DATA  IFIRST/0/ 

SAVE  IFIRST,GAMCUT,GAMMAX 
PARAMETER  (G=9.80D0) 

C 

C***  THE  FOLLOWING  CODE  SEGMENT  COMPUTES  CONSTANTS  REQUIRED  BY  VENTCF. 
C***  IT  IS  EXECUTED  THE  FIRST  TIME  VENTCF  IS  CALLED. 

C 

IF(IFIRST.EQ.O)THEN 
IFIRST  = 1 

GAMCUT  = (2.0D0/(GAM+1.0D0))**(GAM/(GAM-1.0D0)) 

ZZZ=GAM* ( ( 2 . ODO/ (GAM+1 . ODO ) ) ** ( (GAM+1 . ODO) / (GAM- 1 . ODO) ) ) 
GAMMAX  = 3QRT(ZZZ) 

ENDIF 

C 

C***  1.  AND  2.  CALCULATE  THE  P(I)  , DELP,  THE  OTHER  PROPERTIES 
C***  ADJACENT  TO  THE  TWO  SIDES  OF  THE  VENT,  AND  DELDEN. 

C 

DO  10  I = 1,  2 

IF(YREF(I)  . LE.YVENT. AND. YVENT.LE.YLAY(I)) THEN 
C 

C***  THE  VENT  IS  AT  OR  BELOW  THE  REFERENCE  ELEVATION  IN  SPACE  I.  IF 
C***  SPACE  I IS  AN  INSIDE  ROOM  THEN  BOTH  THE  VENT  AND  THE  LAYER 
C***  INTERFACE  ARE  AT  THE  FLOOR  ELEVATION. 

C 

DP(I)  = - G*DENL(I)*(YVENT  - YREF(I)) 

ELSE 

C 

C***  THE  VENT  IS  ABOVE  THE  REFERENCE  ELEVATION  IN  SPACE 
C***  IF  SPACE  I IS  AN  INSIDE  ROOM  THEN  THE  VENT  IS  AT 
C***  CEILING . 

C 

DP(I)  - - G*DENL(I)*(YLAY(I)-YREF(I)) 

$ - G*DENU(I)*(YVENT-YLAY(I)) 

ENDIF 

P(I)  = DPREF(I)  + DP(I)  + PDATUM 
10  CONTINUE 
C 

C***  DELP  IS  PRESSURE  IMMEDIATELY  BELOW  THE  VENT  LESS  PRESSURE 
C***  IMMEDIATELY  ABOVE  THE  VENT. 

C 
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DELP  » (DPREF(2)-DPREF(1))+(DP(2)-DP(1)) 

DO  30  I = 1,  2 

DEN(I)  = O.ODO 
T(I)  = O.ODO 
DO  22  K = 1,  NPROD 

C(K,I)  = O.ODO 
22  CONTINUE 

IF(((YVENT.EQ.YREF(I)) .AND. (YLAY(I) . EQ.YREF(I) ) ) .OR. 

$ ((YVENT.EQ.YCEIL(I)) .AND. (YLAY(I) . LT . YCEIL(I) ) ) )THEN 

DEN(I)  = DENU(I) 

T(I)  = TU(I) 

DO  24  K = 1,  NPROD 

C(K,I)  = CONU(K,I) 

24  CONTINUE 

ELSE 

DEN(I)  = DENL(I) 

T(I)  = TL(I) 

DO  26  K = 1,  NPROD 

C(K,I)  = CONL(K,I) 

26  CONTINUE 

ENDIF 

30  CONTINUE 
C 

C***  DELDEN  IS  DENSITY  IMMEDIATELY  ABOVE  THE  VENT  LESS  DENSITY 
C***  DENSITY  IMMEDIATELY  BELOW  THE  VENT 
C 

DELDEN  = DEN(I)  - DEN(2) 

C 

C***  3.  CALCULATE  VST(I) , THE  "STANDARD"  VOLUME  RATE  OF  FLOW 

C***  THROUGH  THE  VENT  INTO  SPACE  I 
C 

C***  CALCULATE  VST(I)  IF  DELP  = 0 
C 

IF(DELP . EQ . 0 . 0D0)THEN 
VST(I)  = O.ODO 
VST(2)  = O.ODO 

ENDIF 

C 

C***  CALCULATE  VST(I)  FOR  NONZERO  DELP 
C 

IF(DELP . GT . 0 . 0D0)THEN 
VST(2)  = O.ODO 
RHO  = DEN(2) 

EPS  = DELP/P(2) 

ENDIF 

IF(DELP . LT . 0 . 0D0)THEN 
VST(l)  = O.ODO 
RHO  = DEN(l) 

EPS  = -DELP/P(1) 

ENDIF 

X = I.ODO  - EPS 

COEF  = 0.68D0  + 0.17D0*EPS 
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XO  = l.ODO 

EPSCUT  = EPSP*MAX(XO,DPREF(l) ,DPREF(2)) 

EPS CUT  - SQRT( EPSCUT) 

SRDELP  = SQRT(ABS(DELP)) 

FNOISE  = l.ODO 

IF((SRDELP/EPSCUT) . LE . 130 . DO)THEN 

FNCISE  = l.ODO  - EXP (-SRDELP /EPS CUT) 

C 

C***  NOTE:  IF  SINGLE  PRECISION  THEN  USED  65.  INSTEAD  OF  130. 
C 

ENDIF 

IF(EPS.LE.0.1D-  :/)THEN 

W = l.ODO  - 0.75D0*EPS/GAM 

ELSE 

I F ( EPS . LT . GAMCUT ) THEN 

GG  = X**(1.0D0/GAM) 

FF  = SQRT((2.0D0*GAM/(GAM-1.0D0))*GG*GG* 

$ (l.ODO-X/GG)) 

ELSE 

FF  = GAMMAX 

ENDIF 

W = FF/S"^^T(EPS+EPS) 

ENDIF 

V = FNOISE*COEF*W*SQRr(2.0DO/RHO)*AVENT*SRDELP 
IF(DELP.GT.0.0D0)VST(1)  = V 
IF(DELP.LT.0.0D0)VST(2)  = V 
C 

C***  4.  CALCULATE  VEX,  THE  EXCHANGE  VOLUME  RATE  OF  FLOW 

C***  THROUGH  THE  VENT 

C 

C***  STABLE  CONFIGURAT:  M,  CALCULATE  VEX  = 0 
C 

IF(DELDEN.LE.0.0D0)VEX  = O.ODO 
C 

C***  UNSTABLE  CONFIGUPJiTION , CALCULATE  NONZERO  VEX 
C 

IF(DELDEN . GT . 0 . ODO . AND . AVENT . NE . 0 . 0D0)THEN 
IF (NSHAPE . EQ . 1 ) THEN 

CSHAPE  - 0.754D0 

ELSE 

CSHAPE  = 0.942D0 

ENDIF 

DELPFD  - CSHAPE**2*G*DELDEN*D**5/(2.0D0*AVENT**2) 
DPDDPF  - ABS(DELP/DELPFD) 

VEXMAX  - 0.1D0*SQRT(2.0D0*G*DELDEN*SQRT(AVENT**5) 
$ /(DEN(1)+DEN(2))) 

VEX  - MAX (VEXMAX* (l.ODO  - DPDDPF), 0) 

ELSE 

VEX  = O.ODO 

ENDIF 

C 

C***  5.  CALCULATE  VVENT(I)  , THE  VOLUME  RATE  OF  FLOW  THROUGH 


VENTCF  - 25 


c***  THE  VENT  INTO  SPACE  I 

C 

DO  40  I = 1,2 

WENT(I)  = VST(I)  + VEX 
40  CONTINUE 

C 

C***  6.  CALCULATE  THE  VENT  FLOW  PROPERTIES 

C 

DENVNT(l)  - DEN(2) 

DENVNT(2)  = DEN(l) 

TVENT(l)  - T(2) 

TVENT(2)  = T(l) 

DO  50  K - 1,  NPROD 

CVENT(K,1)  = C(K,2) 

CVENT(K,2)  - C(K,1) 

50  CONTINUE 

C 

C***  7.  CALCULATE  THE  VENT  MASS  FLOW  RATES 

C 

DO  60  I - 1,2 

XMVENT(I)  - DENVNT(I)*VVENT(I) 

60  CONTINUE 

C 

C***  8.  CALCULATE  THE  REST  OF  THE  VENT  FLOW  RATES 

C 

DO  70  1=1,2 

QVENT(I)  = XMVENT(I)*CP*TVENT(I) 

DO  65  K = 1,  NPROD 

PVENT(K,I)  = XMVENT(I)*CVENT(K,I) 

65  CONTINUE 

70  CONTINUE 

C 

C***  9.  CALCULATE  THE  RATE  AT  WHICH  THE  VENT  FLOWS  ADD  MASS, 

C***  ENTHALPY,  AND  PRODUCTS  TO  THE  LAYERS  OF  THE  SPACES  FROM 

C***  WHICH  THEY  ARE  EXTRACTED.  FIRST  TREAT  LAYERS  OF  SPACE  1 

C***  (I*=l,  J=2)  AND  TOEN  SPACE  2 (1=2,  J=l)  . 

C 

DO  95  I = 1,  2 

IF  (I.EQ.l)THEN 

J = 2 

ELSE 

J = 1 

ENDIF 

IF ( ( ( ABS ( YVENT - YREF ( I ) ) . LT . 0 . OOOIDO ) 

$ .AND. (ABS(YLAY(I)-YREF(I)) .LT.O. OOOIDO)) .OR. 

$ ( (ABS (YVENT- YCEIL(I)) .LT.O. OOOIDO) .AND. (YLAY(I) 

$ .LT.YCEIL(I))))THEN 

XMU(I)  = -5U1VENT(J) 

XML(I)  = O.ODO 
QU(I)  = -QVENT(J) 

QL(I)  = O.ODO 
DO  75  K = 1,  NPROD 
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PU(K,I)  = -PVENT(K,J) 

PL(K,I)  - O.ODO 
CONTINUE 

IF(ABS (YREF(I) -YCEIL(I) ) .LT, 0 . 0001D0)THEN 
XML(I)  = XMU(I) 

QL(I)  = QU(I) 

DO  80  K = 1,  NPROD 

PL(K,I)  = PU(K,I) 

CONTINUE 

ENDIF 
. ELSE 

XML(I)  = -XMVENT(J) 

XMU(I)  = O.ODO 
QL(I)  = -QVENT(J) 

QU(I)  = O.ODO 
DO  85  K = 1,  NPROD 

PL(K,I)  = -PVENT(K,J) 

PU(K,I)  = O.ODO 
CONTINUE 

IF ( ABS ( YREF ( I ) - YCEIL( I ) ) . LT . 0 . OOOIDO ) THEN 
XMU(I)  = XML(I) 

QU(I)  - QL(I) 

DO  90  K = 1,  NPROD 

PU(K,I)  = PL(K,I) 

CONTINUE 

ENDIF 

ENDIF 

CONTINUE 

RETURN 

END 
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